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Chiral solitons with baryon number one are investigated in the Spectral Quark Model. In this 
model the quark propagator is a superposition of complex-mass propagators weighted with a suitable 
spectral function. This technique is a method of regularizing the effective quark theory in a way 
preserving many desired features crucial in analysis of solitons. We review the model in the vacuum 
sector, stressing the feature of the absence of poles in the quark propagator. We also investigate in 
detail the analytic structure of meson two-point functions. We provide an appropriate prescription 
for constructing valence states in the spectral approach. The valence state in the baryonic soliton is 
identified with a saddle point of the Dirac eigenvalue treated as a function of the spectral mass. Due 
to this feature the valence quarks never become unbound nor dive into the negative spectrum, hence 
providing stable solitons as absolute minima of the action. This is a manifestation of the absence of 
poles in the quark propagator. Self-consistent mean field hedgehog solutions are found numerically 
and some of their properties are determined and compared to previous chiral soliton models. Our 
analysis constitutes an involved example of a treatment of a relativistic complex-mass system. 
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I. INTRODUCTION 

The original proposal by Skyrme of the early sixties 
foresaw that baryons could be described as classical topo- 
logical solitons of a specific non-linear chirally symmet- 
ric Lagrangean in terms of meson fields with non-trivial 
boundary conditions [lj. Within the accepted QCD 
framework the large-iV c analysis under the assumption 
of confinement [EEl supports many aspects of this view. 
Moreover, the identification of the topological winding 
number as the conserved baryon number of quarks from 
the occupied Dirac sea in the background of large and 
spatially extended pion fields Q suggested the under- 
lying necessary fermionic nature of the soliton at least 
for odd N c [5]. The form of the Lagrangean remains 
unspecified besides the requirement of chiral symmetry, 
leaving much freedom on the assumed meson-field dy- 
namics needed for practical computations of low-energy 
baryonic properties d 0] , usually organized in terms of 
a finite number of bosons of increasing mass. One ap- 
pealing feature of the Skyrme model (for reviews see e.g. 
jUSIHilu]]) is that confinement appears to be explicitly 
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incorporated, since the baryon number topology of the 
meson fields cannot be changed with a finite amount of 
energy and the soliton is absolutely stable against decay 
into free quarks. This has also made possible the study of 
excited baryons (see e.g. Ref. [l2| and references therein). 
However, when the partonic content of the baryon is an- 
alyzed in deep inelastic scattering, in the Skyrme model 
the Callan-Gross relation is violated, hence the partons 
turn out not to be spin one- half objects [l3|. Therefore, 
making credible non-perturbative estimates of high en- 
ergy properties is out of reach of the standard Skyrme 
model. This also suggests that the relevant degrees of 
freedom should explicitl y in clude constituent quarks chi- 
rally coupled to mesons [3 EE El El El, OS lH El HI 
motivating the use of chiral quark models with quarks 
and to search for solitonic solutions to describe baryons. 
The approach is very much in the spirit of the Skyrme 
model but with the important feature that the partonic 
interpretation corresponds to spin one-half constituents, 
despite the subtleties of the regularization [23] . The con- 
siderable effort exerted to describe low-lying baryons as 
solitons of effective chiral quark models has been de- 
scribed in detail in the reviews [H, El El El EH • 

The chiral quark models that arise naturally in sev- 
eral approaches to low-energy quark dynamics, such as 
the instanton-liquid model [29( or the Schwinger-Dyson 
resummation of rainbow diagrams (30j . are nonlocal, i.e. 
the quark mass function depends on the quark virtual- 
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ity. For the derivations and applications of these models 

^e^.^MMM H IS Hi M M M IS El, S3, 

EE S, \M, S3]. A major success was the finding 
of baryon hedgehog solitons in nonlocal models [48|, [4t| , 
which turn out to be stable also in the linear version 
of the model. Moreover, the nonlocal models have nice 
features as compared to the local variants, in particu- 
lar they use a uniform regularization of both the normal 
and abnormal parity processes, thus reproducing prop- 
erly the anomalies in the presence of regularization. In 
addition, the theory is made finite to all orders of pertur- 
bation theory. The price to pay for the non-locality is a 
complicated nature of interaction vertices with currents, 
where the gauge invariance imposes non-local contribu- 
tions. Their presence leads to technical complications in 
the treatment of non-local models. 

A relevant issue for chiral quark solitons is related to 
the confinement and stability of the mean field solution. 
The chiral quark soliton model builds baryons as bound 
states of valence constituent quarks in non-trivial me- 
son fields. Indeed, such a model interpolates between 
the non-relativistic quark model and the Skyrme model 
in the limit of small solitons and large solitons, respec- 
tively [50]. As a matter of fact, from the view point of 
the chiral quark model, the apparent confinement in the 
Skyrme model is reinterpreted as a strong binding effect, 
and indeed deeply bound states are not sensitive to the 
confinement properties of the interaction. Nevertheless, 
it is notorious that depending on the details, the mean 
field soliton may in fact decay into three free constituent 
quarks. Moreover, in local models for small enough soli- 
tons the valence state becomes an unbound free quark at 
rest. Actually, to our knowledge there is no satisfactory 
calculation of excited baryon states in the chiral quark 
soliton model precisely because of the lack of confine- 
ment [5l| . These features happen for realistic parameter 
values and despite the alleged compatibility of the chiral 
quark model with the large- N c limit. As already men- 
tioned the soliton description of baryons in the large- N c 
limit is based on the assumption of confinement. 

Let us define the scope of the present work. After 
the many years, the issue of color confinement remains 
a crucial and intriguing subject for which no obvious so- 
lution exists yet. In particular, its realization implies 
two relevant consequences. In the first place, there ex- 
ists a spectrum of color singlet excited states. Secondly, 
quarks cannot become on shell and hence quark propa- 
gators merely cannot have poles on the real axis. This 
restrained meaning of confinement is often called ana- 
lytic confinement in the literature and in that sense is 
adopted in the present work. We hasten to emphasize 
that although our work obviously does not suggest how 
the problem of color confinement might be resolved, we 
manage to set up a framework where the calculation of 
color singlet states such as excited baryons becomes pos- 
sible as a matter of principle within a chiral quark soliton 
approach. 

In this paper we show how a recently proposed ver- 



sion of the chiral quark model, the _Spectral Quark Model 
(SQM) [HI, [H, [Hi], (see also [la] for the original in- 
sights on the spectral problem), not only allows for soli- 
tonic solutions but due to its unconventional and indeed 
remarkable analytic properties yields a valence eigen- 
value which never becomes unbound. Thus the soliton 
is absolutely stable. The model is based on a general- 
ized Lehmann representation where the spectral func- 
tion is generically complex, involving a continuous su- 
perposition of complex masses. The subject of defining 
a well founded quantum field theory of fixed complex 
masses is an old story IH, [H, [6(J HH. Our ap- 
proach differs from these works, since our complex mass 
is an integration variable so that many standard objec- 
tions are sidestepped by choosing an appropriate inte- 
gration contour. In practical terms the spectral function 
acts as a finite regulator fulfilling suitable spectral con- 
ditions but with many desirable properties, in particu- 
lar the simple implementation of chiral symmetry and 
gauge invariance. One of its outstanding features is the 
uniform treatment of normal and abnormal parity pro- 
cesses ensuring both finiteness of the action and a si- 
multaneous implementation of the correctly normalized 
Wess-Zumino-Witten term 1 . This uniform treatment of 
regularization has some impact on soliton calculations 
even in the SU(2) case (where the Wess-Zumino-Wittcn 
term vanishes) since the valence and sea contributions 
to the soliton energy are indeed related to the abnormal 
and normal parity separation of the effective action, re- 
spectively. The basics of SQM are described in detail in 
[531 l55j and the reader is referred there for the descrip- 
tion of the method and numerous applications to the pion 
phenomenology. The SQM approach bares similarities to 
nonlocal models, however the construction of interaction 
vertices with currents is very simple in SQM, as opposed 
to almost prohibitive complications of the nonlocal ap- 
proach. 

The purpose of this paper is twofold: Firstly, we show 
that absolutely stable baryon solitons in SQM exist and 
discuss their properties. Secondly and more generally, we 
study an instance of an involved complex-mass system, 
and show how to treat valence states. Despite the com- 
plexity, the resulting prescription turns out to be very 
simple and easy to implement in practical calculations. 
It amounts to locating the saddles of the valence eigen- 
value, eo, as a function of the complex mass u>, based on 
the condition 

^M = o. (i) 

aw 

The outline of our paper is as follows: In Sect. HT1 we 
provide an operational justification of the need of some 
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part of the Euclidean action and not regularizing the imaginary 
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uniform regularization for the full action without an ex- 
plicit separation between normal and abnormal parity 
contributions. Certain a priori field-theoretic consistency 
conditions are discussed. SQM is shown to fulfill the 
one-body consistency condition in contrast to previous 
local versions of chiral quark models. We also discuss 
in some detail the analytic structure of the meson corre- 
lators, showing that despite analytically-confined quarks 
they possess cuts at large values of q 2 , as requested by 
(asymptotic) quark unitarity. However, the meromor- 
phic structure expected from general large N c consider- 
ations is violated. This point is analyzed in the light of 
a large-iV c Regge model. Next, we pass to the analy- 
sis of solitons. First, we discuss the connection between 
the baryon and topological currents in the limit of large 
solitons, which holds in the presence of the spectral reg- 
ularization. Following the standard field-theoretic ap- 
proach the construction of the baryon state is pursued 
in Sect. IIIII The unconventional appearance of complex- 
masses requires a demanding mathematical treatment of 
both the valence and sea contributions. Nevertheless, 
ready-to-use formulas for the total soliton mass are de- 
rived and analyzed for several soliton profiles showing the 
existence of chiral solitons with baryon number one. An 
alternative derivation is provided in Appendix [E] based 
on computing the spectral integral exactly. Dealing with 
complex mass Dirac operators both for bound states and 
continuum states is involved and some of the features 
may be studied in a somewhat comprehensive toy model 
in Appendix |D] The techniques are close to the more 
familiar complex potentials in non-relativistic quantum 
mechanics which are reviewed for completeness in Ap- 
pendixO In Appendix[G]we show that a linear extension 
of SQM leads to instability of the vacuum, hence SQM 
can only be constructed in the originally proposed non- 
linear realization. In Sect. II VI we look for self-consistent 
hedgehog solutions and determine their properties both 
in the chiral limit as well as for finite pion masses. Fi- 
nally, in Sect. [V] we come to the conclusions. 



II. THE SPECTRAL QUARK MODEL AND 
CONSISTENCY RELATIONS 



We begin with some basic expressions of the general 
field-theoretic treatment of chiral quark models, which 
are the groundwork for our treatment of solitons and the 
method of including valence states in SQM discussed in 
Sect. IIIII In this section we highlight an important con- 
sistency condition which was not fulfilled in the hitherto 
extensively used chiral quark soliton models based on 
local interactions. Remarkably, this condition happens 
to be satisfied in SQM. The analytic properties of the 
quark propagator and in particular the lack of poles are 
reviewed. We also analyze some important aspects of the 
meson sector, and in particular the analytic structure of 
two-point correlators in the complex g 2 -plane. Finally, 
we show how the correct topological current arises in the 



presence of regularization. 



A. A consistency condition 

The vacuum-to-vacuum transition amplitude in the 
presence of external bosonic (s,p,v,a) and fermionic 
(i], fj) fields of a chiral quark model Lagrangian can be 
written in the path-integral form as 



z \j,V,v] = (0|Texp<h 



d A x 



qjq + r]q + qri 



}|0), (2) 



where the compact notation 

j = j> + ^75 - (s + i-f 5 p) 



(3) 



has been introduced. The symbols s, p, v^, and denote 
the external scalar, pseudoscalar, vector, and axial flavor 
sources, respectively, given in terms of the generators of 
the flavor SU (3) group, 



£ ■ 

a=0 



Aa 

1 2 : 



(4) 



with A a representing the Gell-Mann matrices. Any phys- 
ical matrix element can be computed by functional dif- 
ferentiation with respect to the external sources. 

Let us consider the calculation of a bilinear quark op- 
erator, such as, e.g., the quark condensate (for a single 
flavor) (qq). We can think of two possible ways of mak- 
ing such a computation, namely via coupling of an exter- 
nal scalar source s(x) (a mass term) and differentiating 
with respect to s(x), or by calculating the second func- 
tional derivative with respect to the Grassmann external 
sources rj(x) and fj(x) taken at the same spatial point. 
The consistency of the calculation requires the following 
trivial identity for the generating functional: 



(q(x)q(x)) = i 



1 SZ 



Z Ss(x) 



S 2 Z 



lim (q(x')q(x)) = lim (-i) — 
x'^x x'^x Z orj{x)rj{x' ) 



(5) 



where |o means all external sources set to zero. This re- 
quirement can be generalized to any quark bilinear with 
any bosonic quantum numbers and thus we call it the 
one-body consistency condition. In the traditional treat- 
ment of local chiral quark models the l.h.s. of the above 
formula corresponds to a closed quark line and is diver- 
gent, calling for regularization, whereas for x' ^ x the 
r.h.s. is finite and corresponds to an open quark line, thus 
no regularization is demanded. This poses a consistency 
problem which actually becomes crucial in the analysis 
of high-energy processes and introduces an ambiguity in 
the partonic interpretation as well as conflicts with gauge 
invariance and energy-momentum conservation (see, e.g., 
Ref. [12] for a further discussion on these subtle but rel- 
evant issues). Obviously, the previous argument could 



4 



equally be applied to sources with any quantum num- 
bers suggesting that any consistent regularization should 
be applied to the full action. As already mentioned in the 
Introduction, in the traditional approach to local mod- 
els the treatment of singularities requires first to sepa- 
rate the normal and abnormal parity contributions and 
to regularize only the normal parity piece. In the Eu- 
clidean space this separation corresponds to the real and 
imaginary parts of the action. 

In the next section we show that SQM fulfills the 
one-body consistency condition. More generally one 
might want to extend the condition ([5]) to any number of 
quark bilinears contracted to bosonic quantum numbers, 
namely the N— body consistency condition 

{q{x-i)Txq{xi) ■ ■ ■ q(x N )I \q(x N )) 
= lim (q(x' 1 )T 1 q(x 1 )---q(x' N )Txq(x N )) (6) 

x'. — >Xi 

where the l.h.s. is evaluated after functional derivatives 
with respect to bosonic sources and the r.h.s. with re- 
spect to fermionic sources respectively. Here Tj are gen- 
eral spin-flavor matrices. We have found that this is 
not possible in the SQM scheme (see the discussion in 
Sect. HH below). 



B. The Spectral Quark Model and the quark 
propagator 

In SQM the regularization is imposed already at the 
level of one open line through the use of a generalized 
Lehmann representation for the full quark propagator in 
the absence of external sources, 



and 



duj- 
c V 



(7) 



where C is a suitable contour in the complex mass- 
plane (see also Fig. 1 below). Chiral and electromag- 
netic gauge invariance can be taken care of through the 
use of the gauge technique of Delbourgo and West [63|, 
[64]], which provide particular solutions to chiral Ward- 
Takahashi identities, or through the use the standard 
effective action approach, applied in this paper. As a 
result, one can "open the quark line" from one closed 
loop and compute high-energy processes with a partonic 
interpretation, such as, e.g., the structure function and 
the light-cone wave function of the pion, etc. 15411 . or the 
photon and p- meson light cone wave functions [651 ] . The 
generalization of the (one body) consistency condition 
for scalar sources, Eq. |(SJ), to all possible bosonic quan- 
tum numbers can be achieved by taking the generating 
functional to be 



Z[r),r),s,p, . 



DJJ e -i{v,S[U,s,p,v,a]r)) e iT[U,s,p,v,a] ^ 



where the quark propagator and the effective action are 
given by 

(x'\S[U )S) p,v,aU>\x)= f dup(uj)(x\(D)~^\x'), (9) 

Jc 



T[U, s,p, v, a] = ~iN c / dwp(u>)Tr log (iD) , (10) 
Jc 

respectively. The Dirac operator has the form 

iD = i$ — ujU 5 — too + + ^75 — s — ij 5 p) ■ 

For a bilocal (Dirac- and flavor-matrix valued) operator 
A(x,x') we use the notation 

TtA = [ d i xtr(A{x,x)), (11) 



the flavor trace. 



with tr denoting the Dirac trace and 
The matrix 

= \{l +15P + 1(1-75)^, (12) 

while U — u 2 — e lv/ ^*/^ is the flavor matrix representing 
the pseudoscalar octet of mesons in the nonlinear repre- 
sentation, 







7T+ 




K+ 


7T 










K~ 




K° 







(13) 



The matrix too = diag(TO u , to^, to s ) is the current quark 
mass matrix and / = 86 MeV denotes the pion weak- 
decay constant in the chiral limit, ensuring the proper 
normalization condition of the pseudoscalar fields. For a 
two-flavor model it is enough to consider $ = f ■ if/ a/2- 
In Eq. ([8]) the Dirac operator appears both in the 
fermion determinant as well as in the quark propagator 
and could, in principle, be treated independently. The 
one-body consistency condition is fulfilled precisely be- 
cause we use the same spectral function p(u>) for both. 
We have refrained on purpose from writing a Lagrangean 
in terms of quarks explicitly since anyhow chiral quark 
models are defined in conjunction with the regularization 
and the approximation used. In our case we work in the 
leading order of the large- N c expansion, which amounts 
to a saddle point approximation in the bosonic [/-fields 
and use the spectral regularization which is most explic- 
itly displayed in terms of the generating functional pre- 
sented above. The new ingredient of SQM compared to 
earlier chiral quark models is the presence of the quark 
spectral function p(uj) in Eq. (|9I10[) and the integration 
over u> along a suitably chosen contour C in the complex- 
ly plane. Our approach extends the early model of Efi- 
mov and Ivanov [56j by including the gauge invariance, 
the chiral symmetry, and the vector meson dominance, as 
well as applying the model to both low- and high-energy 
processes. 

The above SQM construction implements the one-body 
consistency condition, as follows 



SZ 



= lim (— i)' 



s 



5r)(x) r)(x') 



(14) 
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which is more general than Eq. |(5j since it is valid in the 
presence of non-vanishing external bosonic sources with 
any quantum numbers and non-trivial background pion 
field, U . This represents a clear improvement on the pre- 
vious local chiral quark models where this requirement 
was violated. Eq. (fT4|) has direct applicability on the 
soliton sector as we will discuss in Sect. IIII1 It should be 
mentioned, however, that similarly to many other chiral 
quark models, the two-body and higher consistency con- 
ditions, Eq. © are not satisfied in SQM. For instance, 
for the two body consistency condition one gets 

(q(xi)Tiq(xi)q(x 2 )T 2 q(x 2 )) = — J DU J dujp(uj) 



x tr [r^nlD-V)!^}] tr [r 2 (z 2 |D-V)|z2) 



while 



lim (q(xi)Tiq(x' 1 )q(x2)T 2 q{x' 2 )) = 

X . — >Xi 

x / dwip(wi)tr [r^ailD-^wiJIari)] 



dw 2 p(uj 2 )tr [ra^lD" 1 ^)^!)] 



^ {q{xi)T 1 q{x 1 )q{x 2 )T 2 q{x 2 )) 



lim - 



(15) 



DU 



(16) 



Likewise, we do not know if this one body consistency 
condition is fulfilled beyond the leading large- N c approx- 
imation. At present the only known way to fulfill all 
consistency conditions is by returning to nonlocal ver- 
sions of the chiral quark model. Therefore, we must pro- 
vide a prescription of how higher functional derivatives 
should be handled in case were more than a single quark 
line could be opened. Eq. (fl4|) suggests to use always 
the method based on the bosonic sources for operators 
involving any number of quark bilinear and local opera- 
tors. For operators involving one single bilocal and bi- 
linear operator q(x)q(x') we may use the method based 
on fermionic sources, since Eq. (fTi)) guarantees the con- 
sistency in the coincidence limit x' — > x. Of course, this 
prescription does not yield a unique result when more 
then one bilocal and bilinear quark operator is involved 
or equivalently when more than one quark line is opened. 

Using the standard variational differentiation for the 
generating functional the Feynman rules for SQM follow 
from the action (flO]) . They have the form of the usual 
Feynman rules for a local theory, amended with the spec- 
tral integration associated to each quark line, according 
to Eq. (|10[) . This resembles very much the well known 
Pauli-Villars regularization method (however with a con- 
tinuous superposition of complex masses) and allows for 
very efficient computations, see (52l. l53l. [54. l55l| . 

The basic paper [53| explains the general construction 
and the conditions for the moments of the spectral func- 
tion p(u>) coming from physics constraints. In particular, 
normalization requires 



dojp(uj) = 1, 



(17) 



while observables are related to the log-moments and in- 
verse moments of p(u>) 53]. The full spectral function 
consists of two parts of different parity under the change 
lu — > — u), i.e. p(u>) — pv(to) + Ps{u), with the scalar 
part, ps, even and the vector part, py, odd. A par- 
ticular implementation of SQM is the meson dominance 
model, where one requests that the large- N c pion elec- 
tromagnetic form factor has the monopole form of the 
vector-meson dominance (VMD), 



FAq 2 ) 



m 2 



(18) 



where My denotes the p-meson mass. 2 The matching 
of the model predictions to this form yields the rather 
unusual spectral function 



Pv{uj) 



1 1 



1 



2niuj (1- Auj 2 /M 2)5/2 = 



(19) 



with the pole at the origin and cuts starting at ±My/2, 
where My is the mass of the vector meson. Similar con- 
siderations for the photon light-cone wave function [ssj ] 
and matching to VMD yield for the scalar part an anal- 
ogous form, 



1 



48tt 2 



99 



2m N C M S (1 ~Auo 2 /M 2 f/ 2 ' 



(20) 



where (qq) is the single flavor quark condensate, and 
A/5 = My (55[. The contour C for the to integration 
to be used with formulas (| 1 9I20|) is shown in Fig. Q] for 
pv (for ps one has the same contour with My — > Ms). 
This contour is applicable in the vacuum (i.e., no baryon 
number) sector of the model. The extension to baryons 
is described in the next section. From Eq. (fl"9|) and ([20]) 
one gets the quark propagator functions from Eq. ([7]), 



A(p 2 



B(p 2 



p v (oj) 
duj— 

c P - w 



1 



1 

dio 



(1- V/M 2 ) 5 / 2 
u>p s {u>) 



c P 2 - 



48ir 2 (qq) 



M S N C {1-Ap 2 /M 2 f/ 2 ' 



(21) 



These functions do not possess poles (the alleged pole at 
p 2 = in A(p 2 ) is canceled), but only cuts starting at 



2 This example shows in a transparent way a peculiar feature of 
the model. In the standard constant mass case the pion form 
factor, a three point function, due to Cutkosky's rules, displays 
a cut in the form factor for sufficiently large energy due to a 
superposition of poles in the quark propagator. In SQM the 
mechanism is just the opposite; the cuts conspire to build a pole 
in the form factor. 
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-Mv/2 



w-Complex Plane 



Mv/2 



FIG. 1: The contour C in the complex-o; plane for the calcu- 
lations in the vacuum sector in the meson dominance variant 
of SQM. My denotes the p-meson mass. The two segments 
shown in the figure are connected at ± infinity with semi- 
circles, not displayed. 



p 2 = My/4, reflecting the structure of p(ui). Actually, 
the cj-integral can be evaluated using the integral trans- 
formations displayed in the first lines of Eqs. pTjl . For 
instance, taking the limit p 2 — > oo one gets the moments 



du py(oj)u) 



2k 



fc = 0,l,2,. 



(22) 



c 



= / duo Ps {uj)uJ 2k+1 , k = 0,1,2,... (23) 
Jc 

which imply in particular the vanishing of all positive 
moments for the spectral function p(u>). 

The pion weak-decay constant in the chiral limit comes 
out to be [531 



2 _ N C M 2 
1 24tt 2 • 



(24) 



a relation which works well phenomenologically. This 
relation will be used as an identity in the rest of the 
paper. Compared to the standard field-theoretic case, 
each quark line is supplemented with a spectral integra- 
tion f c du>p(uj). This makes calculations very straight- 
forward and practical for numerous hadronic processes 
involving mesons and photons. We stress that despite a 
rather "exotic" appearance of the quark spectral func- 
tion, SQM leads to proper phenomenology for the soft 
pion, including the Gasser-Leutwyler coefficients [55[ , the 
soft matrix elements for hard processes, such as the dis- 
tribution amplitude, transition form factor, or structure 
functions [53j, t he g eneralized forward parton distribu- 
tion of the pion 67| , the photon distribution amplitude 
and light-cone wave functions 16511. o r the pion-photon 
transition distribution amplitude [661 ] - In addition, the 
calculations are straightforward, leading to simple an- 
alytic results. Interestingly, the quark propagator cor- 
responding to (|19l20p has no poles, only cuts, in the 
momentum space. The evaluated mass function of the 
quark, M{Q 2 ), displays a typical dependence on the vir- 
tually Q 2 in the Euclidean region and at the qualitative 



and quantitative level compares favorably to the non- 
local quark models and to the lattice calculations [68[. 
As a sample calculation, we extend in the next subsec- 
tions the previous considerations to the calculation of 
vacuum properties and two-point correlators. 



C. The vacuum sector 

As mentioned above the functional (one-body) consis- 
tency conditions guarantee the unambiguous calculation 
of observables obtained from quark fields, either as dou- 
ble Grassmann functional derivatives or as single bosonic 
ones. These identities are generically formally satisfied 
but become tricky under regularization. Here we use the 
first possibility for the quark condensate and the vacuum 
energy density. The second method based on the effec- 
tive action and yielding identical results is outlined in 
Appendix [A] 

The (single flavor) quark condensate (qq) can directly 
be computed from the quark propagator 



-%N C 



d 4 P 



■TrSV) 



N f (qq) = -u, cJ 

= J <0r* B{p% (25) 



which through the use of Eq. (j2"Tj) becomes an identity. 
The vacuum energy density can be computed by using 
the trace of the energy momentum tensor, yielding 



(?(*)- 



-iN r 



d 4 P 

(2tt)< 



Tr 



I(^V+pV)^) 



= -UN, 
1 



d 4 P 



lNcNf J Wr^ p ^ A(p ) 



(26) 



where the subscript refers to the trivial vacuum (A = 1 
and B = 0). The bag constant is defined by the differ- 
ence e = \{0£) — |(#u)o an d after computing the integral 
becomes 



M 4 N f N c 
192vr 2 



(27) 



The formula implies that it costs energy to dig a hole in 
the vacuum with a non-vanishing pion weak decay con- 
stant (see Eq. (HU)), as one might expect for the chirally 
symmetric broken phase. 



D. Two-point correlators in the meson sector 



In SQM the quark propagator has no poles but cuts 
at p 2 = My/4 preventing the occurrence of quarks on 
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the mass shell. This complies to the notion of analytic 
confinement, a necessary but certainly not sufficient con- 
dition for color confinement. This non-standard behavior 
already suggests a possible departure from the standard 
treatment. However, to justify the claim of confinement 
one should check, in addition, the absence of cuts in phys- 
ical correlators. Note that this question is in principle 
unrelated to the analytic confinement of bound quarks 
in a soliton with baryon number one which will be dis- 
cussed in further sections below, but it is still interesting 
to see the analytic properties of the model. We analyze 
this issue below in some detail for the two-point mesonic 
correlators defined as 



IW<?) 



\T{J A (x)J B (0)}\0), (28) 



where Ja{x) and Jb(x) are interpolating currents with 
the relevant meson quantum numbers. In the standard 
path integral approach (see Eq. {8])) the time-ordered 
products of currents in the vacuum can be evaluated 
by suitable functional derivatives of the generating func- 
tional (JTDJ) with respect to external bosonic currents, 
and in the large- N c limit the path integral is driven to 
the saddle point of the path integral in the presence of 
those currents. On the other hand, at large N c any two- 
point mesonic correlation function should have the gen- 
eral structure [2j, |3| 



n,/(g 2 ) = £ 



(29) 



due to confinement. This is a very stringent test, since 
it implies in particular a meromorphic analytic structure 
in the q 2 complex plane. Amazingly, this rather simple 
requirement has never been accomplished in chiral quark 
models, either local or nonlocal 3 . 

A straightforward consequence of the large- N c repre- 
sentation of the two point correlator in Eq. (|29[) is posi- 
tivity, Hj(q 2 ) > 0, in the Euclidean region q 2 = —Q 2 < 
and quark unitarity at large values of q 2 . This ensures 
that, for instance, the inclusive e + e~ — ► hadrons total 
cross section is proportional to the imaginary part of the 
polarization operator in the vector channel. 4 

In our previous work [53j we have evaluated the VV 
and AA correlators corresponding to the conserved vector 
and axial currents 



A„ 



JA a ( x ) = q( x )Y l i5^-q(x) 



(30) 
(31) 



3 Instead, calculations at large values of Q 2 have only been con- 
fronted to QCD sum rules (see, e.g., Ref. [47l and references 
therein.) 

4 This asymptotic quark unitarity holds for external currents and 
should not be confused with the S-matrix hadron unitarity which 
is violated at any finite order of the large N c expansion due to 
the l/v^Vc suppressed behavior of the meson- meson couplings. 



for a general spectral function p{oj) using solutions to the 
Ward-Takahashi identities based on the gauge technique. 
In this paper we use the effective action (fT0| to evaluate 
the correlation functions. Details of the calculation are 
provided in Appendix [A"l We find 



K a v\<l) = \Sa b (g^q 2 - f<t) ^(q 2 ), (32) 



ith 5 



2N r 



IW) = ^ / p(u)du X 



(33) 



and 



with 



2u 2 (I(q 2 ,u) - 1(0, w)) - q 2 (^ - I(q } . 



KTHq) = (^V " W), (34) 



rr^) = n^(^) + m c J duo^pi^iiq 2 ^). (35) 

In the particular case of the meson dominance model 
for the spectral function, Eq. (fl9|) . up to an (infinite) 
constant one has a remarkably simple result 



/i 



M 2 



- q 2 24tt 2 
2 _N, 

2 ~ 



^-24^ l0g ( v ^ 



(36) 



where fy = f 2 . The above formulas clearly display the 
p meson pole in the vector channel and the pion pole in 
the axial channel. Also, we see that f\ = 0, i.e., there 
is no axial meson dominance. Our expressions fulfill the 
first Weinberg sum rule, 



hm/ [nW)-^W)]=.f 2 



(37) 



Similarly as in other local quark models, the second 
Weinberg sum rule is not satisfied, 

km q 4 [nl(q 2 )-Il T A (q 2 )] =-M 2 f 2 ^0, (38) 

q z — > — oo 

as noted in our previous work [53| . Since the sec- 
ond Weinberg sum rule is a high-energy feature of the 
theory, one hopes that it is not crucial for the low- 
energy phenomena studied in this paper and in particu- 
lar for the soliton properties which probe Euclidean mo- 
menta corresponding to a soliton size ~ y/E/My (see, 
e.g., Eq. (|98[) below). This assumption has implicitly 



We correct here a typo in our previous work [53j |. Eq. (4.3). Our 
conventions here are Tl v (q 2 ) = —2Tl(q 2 )/q 2 . 
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been made in state-of-the-art chiral quark soliton mod- 
els l2l|2ll2&l27L|2l. 



The presence of the log pieces in the correlators guar- 
antees the fulfillment of quark unitarity since 6 



infinite Regge sum 



W(l + z) - *(1) 



£ 

n=l 



1 



n + z 



(40) 



ImH^ 2 ) = fmlT^ 2 ) = 



q 2 > My. 



(39) 



The appearance of these quark unitarity cuts can be in- 
ferred from the general structure of the quark propagator 
(see Eq. (f2Tjl ) for any quark momentum p 2 > My/4, in 
spite of the fact that there are no poles in the quark 
propagator. A more detailed analysis is presented in Ap- 
pendix [Bj Note that in QCD one obtains these parton- 
like relations for q 2 — > 00. We stress that the coefficients 
of the log terms are precisely such as in the one-loop 
QCD calculation, complying to the parton-hadron du- 
ality. Moreover, despite the unconventional features of 
SQM involving complex masses, positivity is preserved 
both for the pole and for the cut contributions to the 
imaginary parts of the considered correlators, as can be 
seen from Eq. (f36|) . This is a virtue of the spectral model, 
not easily fulfilled in other chiral quark models 7 , ft is 
important to realize that this quark unitarity relation 
is hidden in the large- N c meromorphic representation of 
Eq. (l29l) a t large q 2 through the asymptotic density of qq 
states [70j . 

To see this we may compare to the large- iV c Regge 



models (see e.g. [7l| and references therein), where the 



meson spectrum is chosen to be a tower of infinitely many 
radially excited states with masses Ml v = My + 27rc7n, 
where a is the string tension and the residue is taken to 
be constant / 2 v = N c a/(12w) precisely to implement 
quark- hadron duality at large q 2 . As we see, SQM cor- 
responds to keeping one pole in the vector channel and 
approximating the higher states by a logarithm. More 
explicitly, one has 



IIy(g' 

N c a 

71=0 



iV(o) 



t 



My 



24tt 2 



2iran — q 2 

M v -q 2 \ (M* 
27TCT / \2ira 



My + 2nan 



where — ^(1) = 7 is the Euler-Mascheroni constant, by 
the approximation of taking explicitly the first pole. The 
higher poles lead to the asymptotic behavior 



l/(z + l) + l + log(l + z), 



(41) 



approximated by the cut in SQM. The accuracy is bet- 
ter than 20 % for < z < 00. We also recall that the 
large- N c analyses restricted to a finite number of reso- 
nances [72I [73f provide a meromorphic structure but fail 
to give the large-q 2 parton-hadron duality conditions. 

For completeness we also discuss the scalar and pseu- 
doscalar channels, where the currents are given by 



J|(x) = g(i)yg(i), 
Jp(x) = q{x)i-i b ^-q{x). 



(42) 
(43) 



No Ward-Takahashi identities may be written for these 
currents and thus they are not directly amenable to the 
gauge technique. The effective action approach used here 
is superior to the gauge technique since it also allows to 
compute the two-point correlators not directly related to 
conserved currents. The result is given up to two sub- 
traction terms as follows (see Appendix IA"]). Defining 



n a P b P (q) = ^s ab n P ( q ) 



(44) 



one gets 



N, 



2 <? 2 lo S 



N, 

2(qq) 2 



M 2 



My 



V 2 M 2 
My — q 



1 M|(M y -g 2 ) 

P q 2 M v {M 2 -q>y 



(45) 



where the digamma function, fy'(z) — Y'{z)/T{z) has 
been introduced. SQM corresponds to representing the 



We define the discontinuity as 

Discl% 2 ) = n(g 2 + i0+) - U(q 2 - i0+) = 2ilmn( 5 2 ). 



where the residue of the pion pole is 2/ 2 i? 2 , with Bq = 
— (qq)/f 2 - We see that fs — 0, i.e., there is no excited 
pseudoscalar meson dominance. Again, the emergence 
of the quark unitarity cuts is manifest. Similarly to the 
vector and axial correlators, the first SS — PP Weinberg- 
like sum rule is verified 

lim q 2 [n s (q 2 ) - U P (q 2 )] = ~2B 2 f, (46) 



7 For instance, the Pauli-Villars regularization spoils positivity 
due to subtractions, while dispersion relations are fulfilled. The 
proper-time regularization preserves positivity but does not ful- 
fill dispersion relations due to a plethora of complex poles |69H . 



while the second SS — PP sum rule is not, 

lim q 2 [Il s (q 2 ) - H P (q 2 )] = -2M 2 f 2 ? 0. (47) 
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Our analysis of both the second Weinberg sum rules for 
VV — AA and SS — PP correlators agrees with the ob- 
served mismatch of the low-energy coefficients Lg and 
L\o between SQM [55[ and the large- N c evaluation in the 
single-resonance approximation (SRA) [j^, [73| (where 
both sum rules are enforced). In [55| it was also shown 
that matching L3 of both SQM and the large- N c SRA 
yields the identity between the scalar and meson masses, 
A/5 = My (as we assume for the rest of the paper). It is 
interesting that the same condition also follows from the 
requirement of having a single pole in the PP correlator, 
Eq. Moreover, at small q 2 one has 



n v (g)-n A (g) 



r- 



AL 



n s (g)-n P (g) = -2B 2 
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+ 16L 8 

q 2 



(48) 
(49) 



Using Eqs. ([Ml) and (@SJ) and taking M s = M v yields 



N r . 



384tt 2 

N n 



\Q{qq? 



92tt 2 



(50) 



in agreement with results from the derivative expansion 
carried out in 55]. The large- N c SRA yields L s — 
3/ 2 /(32Mj) and L w = -3/ 2 /(8M£) [zllzal, whil e the 
large- N c Regge models produce L w = — N c / (96y/3n). 8 

To summarize this section, SQM provides meson two- 
point correlators which carry poles as well as cuts, while 
strictly speaking large N c requires only having poles. 
Thus, despite the quark propagator not having poles, the 
meson correlations do not exhibit true confinement. Al- 
though we cannot prove it in general for any correlation 
function, all two-point correlators we have considered 
obey positivity and analyticity, i.e. dispersion relations. 
This is required by the (asymptotic) quark unitarity and 
causality, i.e. the current commutators must vanish out- 
side the causal cone, [j(x) 7 j(0)] — for x 2 < 0. The 
coefficients of the log terms are in agreement with the 
parton-hadron duality. On the other hand, the second 
Weinberg sum rule is not satisfied, as in other local chi- 
ral quark models. This calls for caution, in particular 
when analyzing processes sensitive to high-momenta. 



8 At the mean field level in a gradient expansion the Lg term cor- 
responds to C(m 4 ) corrections to the soliton energy, while L10 
couples to external axial and vector currents. So a slight mis- 
match in these values should not influence strongly the mean field 
soliton properties. We estimate the corresponding correction to 
the energy as A_B 8 = AL a 4m* J d 3 x[cos(26»(r)) - 1] ~ -3MeV, 
with ALg ~ 10~ 3 — a negligible number. 

9 We remind that nonlocal models do indeed fulfill this high energy 
constraint [39j and the violation of the second Weinberg sum rule 
in local models is probably related to the violation of the two 
body consistency conditions discussed above (see Sects. Hi Al and 
Sect. nTE) , On the other hand it is uncertain if nonlocal models 
do indeed satisfy analyticity. 



E. The topological current 

In a previous work [55| it was shown how the Wess- 
Zumino-Witten term arises for SQM in the presence of 
the spectral regularization. As a sample calculation illus- 
trating the consistency condition of Sect. Ill Al we com- 
pute the baryon current in the limit of spatially large 
backgrounds. This also shows how the regularization ef- 
fects cancel in the final result, yielding the well known 
Goldstone-Wilczek current Taking the appropriate 
functional derivative in Eq. ([2]) we find 



(qixh^qix)} 



1 SZ 



Z 5v^{x) 
lim (— i) 2 





1 S 



SZ 



r 



Z Srj(x) Sn(x') 



dLop{uj)tv 



-f{x\- 



LOU 5 ' 



.(51) 



The first two lines display the consistency condition. The 
derivative expansion of the Dirac operator can be neatly 
done with the help of the identity 



1 



1 i@ - ujU 5 



\x) = tr 



d 4 k 



1 



(2tt) 4 ft + i@ - ojU 5 ' 



(52) 



where the differential operator acts to the right on the 
function f(x) = 1. This formula can be justified through 
the use of an asymmetric version of the Wigner trans- 
formation presented in Ref. [74| . Formal expansion in 
powers of d yields 



00 



d 4 k 



1 



-1 rc+l 



A' 2 



ujU' ' J (2tt) 4 

n=0 J v ' 

Qt + loU 5 ^) [(-#) (jg + uoU^)] n . (53) 
The leading non- vanishing term is 



(q(x)^q{x)} 



c dUJP{u) J (2^ 
4 



[k 2 



u> 



214 



tr 



(54) 



where -B M (x) is the topological Goldstone-Wilczek cur- 
rent, 



B„ 



1 



24tt 2 



^"((rfdcJJ) (U^dpU) (rfduU)), (55) 



and 



I r>r, 2 \ f r \ f d A k LO A 

Po = (-96. ,) Jj^) J (2?r)4 [fc2 _ ^ 2]4 



duip(uj) = 1, 



(56) 



where in the first line we have computed the momentum 
integral and in the second line used the normalization 
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condition (JTTJ) . Thus the proper normalization of the 
spectral function is responsible for the correct topological 
properties in SQM, preservation of anomalies, etc. In 
fact, we view the uniform treatment of the anomalous 
processes as one of the main advantages of SQM over 
other local chiral quark models. 

III. BUILDING THE BARYON 

In this section we show how valence orbits are con- 
structed in SQM. The issue is far from trivial, as un- 
like the case of local models, there is no obvious Fock 
space representation. Already the experience of nonlo- 
cal models showed that the construction of valenceness 
is an involved issue [HI, H^] . There a bound-state pole of 
the propagator was found in the background of hedgehog 
chiral fields, and this state was occupied with iV c = 3 
valence quarks. The full contribution to the baryon cur- 
rent (local and nonlocal) yielded a correct (and quan- 
tized) baryon number of the soliton. In the present case 
we face a different situation, with the spectral density 
and inherent complex masses present. Thus we start our 
derivation from very basic field-theoretic foundations, ar- 
riving in the end at a very simple prescription holding 
under plausible assumptions concerning the analyticity 
properties of the valence eigenvalue as a function of the 
complex mass. 

A. Interpolating baryon fields 

In the standard field-theoretic approach a baryon can 
be described in terms of the corresponding correlation 
function 

tt B (x, x) = (0\t{b(x)B(x')}\0) (57) 

with B(x) denoting an interpolating baryonic operator 
in terms of anti-commuting quark fields. We take the 
simplest combination 

B{X) = -^e a ^-- a ^^-^q aiai (x) ■ ■ ■ q aNc a Nc (x) , 

(58) 

where (ax, . . . ,o>n c ) are the color indices, (a\, . . . , ajv c ) 
the spinor-flavor indices, and <j> a i> - > a N c j s the appropriate 
completely symmetric spinor-flavor amplitude. Inserting 
a complete set of baryon eigenstates gives 

U B (x,x')= (59) 
6(t - Y,{0\B(0)\B n , k)(B n , k\B(0)\0)e-^ x ^ k 

n 

+(-l) N °e(t' - t) ^(0|B(0)|S„, k)(B n , k\B(0)\0) 

a 

xe +i( x ~ x ') k 



where B n (B n ) are the baryon (antibaryon) states with 
momentum k. Next we take the limit t — t' = T — * — ioo. 
That way the lightest baryon at rest is selected in the 
sum, 

IL B (x,x') = (Q\B(Q)\B)(B\B{0)\0)e- l ™ B . (60) 

In the large- N c limit, we first rewrite the time ordered 
product as a path integral over the fermionic degrees 
of freedom with the weight exp(iSsQM)- The resulting 
expression, in turn, can be obtained by the appropri- 
ate functional differentiation of the generating functional 
Z[s,p, v, a, ?y, fj] with respect to the external quark fields 
r)(x) and fj(x), yielding 

U B (x,x') = $«i,-.«iV c $o' 1 ,...,o' Wc 
f N c 

x DUe ir ^l[iS ai< (x,x'). (61) 
•* i=i 

Here the one-particle Green function in SQM is given by 

S aa ,(x,x')= f (kjp(cj){x\(i^)-ijU 5 )~^\x') (62) 
Jc 

and the one loop effective action has the form 

T[U] = -iN c [ dup{oj)Trlog (i@ - ujU 5 ) . (63) 
Jc 

Note the presence of the spectral integration, however, 
the contour C to be used in the baryon sector is yet to 
be determined. There is no a priori reason why it should 
be equal to the vacuum contour C 10 . The limit N c — * 00 
drives the functional integral over the bosonic U fields 
into a saddle point. In the vacuum sector one then ob- 
tains the mean-field vacuum. In the baryon sector, due 
to the presence of N c factors S aya > in the numerator, the 
dominating saddle point configuration is of course differ- 
ent from that of the vacuum and depends non-trivially 
on x and x' . The limit of a large evolution time T selects 
the minimum-energy stationary configuration. For sta- 
tionary configurations the Dirac operator can be written 
as 

i$-uU 5 =yo(id t -H(uj)), (64) 

with the Dirac Hamiltonian equal to 

H{uj) = -ia- V + luU 5 . (65) 

Note that the spectral mass to appearing here is com- 
plex and the situation is unconventional, since H(u>) is 



Note by analogy that in the standard many-body quantum field 
theory the contour in the complex energy variable selects the 
orbits to be occupied and is clearly different for states with dif- 
ferent baryon number, or baryon density, where it crosses the 
real axis at the value of the chemical potential. 
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not Hermitean. In such a case one must distinguish be- 
tween the right and left eigenvectors, Hijj, 



R _ 



„V„ and 



c 

tri^ni n °t related by the Hermitean conjuga- 



tion, i.e. (ip^) ^ ^n- The orthogonality relation is 
(ipn>' t Pm) = Snm and completeness is given in terms of 
the left-right identity J2 IV'nKV'nl = 1- We will not use 
this fancy notation, but will implicitly understand that ip 
is the right eigenvector, while ip* is in fact the complex- 
conjugated left eigenvector. The corresponding eigen- 
value problem is then 



H(u)<ip n (x;uj) = e n (aj)ip n (x;uj). 



(66) 



For our particular Hamiltonian one has the following 
useful properties 

J 5 H(uj) 7 ^ = H(-u), 

(7o75)^(w)(7o75) _1 = -H(oj), 

trH(u) = 0. (67) 

where the trace is in the Dirac sense. Some properties of 
the eigenvalues deduced from the properties (|67p are 



e n (u>)* = e n (uj*), 
e„(w) = e n (-w). 



(68) 



One can now use the spectral representation of the prop- 
agator 



(69) 



iS aa ' (x, t; x , t) = / dup(u;) 

Jc 

all/ e iu(t-t') 

y 7. 7. A 1pna(x;u)lpna>(x';Uj) 



v - e n (u) 



with ip na (x;Lj) (ip na {x; Co)) and e n (u>) denoting the right 
(left conjugated) eigenfunctions and eigenvalues of H(u) 
evaluated at the stationary bosonic configuration. The 
contour for the energy integration, discussed in the fol- 
lowing, is denoted by 7. For the calculation of the baryon 
state note that the propagator given by Eq. (f6"9"| and the 
fcrmion determinant given by Eq. (|63[) have to be evalu- 
ated at large Euclidean times. 



B. Hedgehog ansatz 

The stationary solutions of chiral quark models have 
the familiar hedgehog form, 



U 5 (x) = exp(i7 5 f • 0(f)), 



(70) 



with denoting the chiral phase field. In the hedgehog 
ansatz 



(71) 
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FIG. 2: The lowest J p = 0+ eigenvalue eo as a function of 
mass M for several soliton size parameters R in the profile 
(|72[) . The curves are arranged with the lowest value of R at 
the top. The two straight solid lines indicate the boundaries 
of the gap, M and —M. Each curve leaves the positive con- 
tinuum at some critical value of M, assumes a maximum, and 
then decreases, asymptotically becoming parallel to the — M 
line. 



with 0(r) being a radial function. In the following two 
subsections we analyze the Dirac problem (|66j) for the 
exponential profile as an example, 



0(r) 



7re 



-r/R 



(72) 



In particular our figures are obtained for that case. Other 
popular cases include the linear profile 



0(r) = w{l-r/R)6(R-r) 
and the arctan profile [2(| 

9(r) = 2arctan(i?/r) 2 . 



(73) 



(74) 



The parameter R is a generic size scale of the soliton 
proportional to the baryon rms radius. The profile (|72p 
is in fact a fairly good approximation to the fully self- 
consistent profile found numerically in Sect. HVl 

In the treatment of hedgehog systems it is relevant to 
consider symmetries such as the grand spin G = I + J, 
the sum of isospin and spin. The reader is referred to 
Refs. [H [H Ha H3] for necessary details. 



C. The valence contribution 

Now we resort to the numerical calculation of the spec- 
trum in the chiral soliton, described in detail in the pro- 
ceeding sections. What we need for the present discussion 
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is the behavior of the spectrum when we vary the spec- 
tral mass u> as an independent variable. We focus on the 
grand-spin parity + state, G p = + , where G = I + J 
is the sum of isospin and spin, appropriate to classify 
states in the hedgehog background (fTTj) . In Fig. [2] we 
show the dependence of the lowest + eigenvalue on the 
mass parameter M for real masses (in our notation AI 
is the same as real to). Subsequent curves are obtained 
for the profile ([72]) with a fixed radius R. In fact, all 
the displayed curves are related to each other via simple 
scaling. Indeed, introducing £ = r/R and £1 = ujR we 
can rewrite the Dirac equation in the form 



[~ia • V f - + nu"(0) fi) = £„(fi)*„(£;ft), (75) 

which depends on a single scale Q. By identification 

E n {uuR) 



R 



(76) 



As can be promptly seen from Fig. [2j for any value of 
R we have for the valence + level the limiting behavior 



e (M) 
eo(M) 



M 



-M + a 



at 
at 



M 
M 



0, 

oo, 



(77) 
(78) 



where a is a constant. Correspondingly, at fixed M in the 
limit of low R (small solitons) we have cq(M) ~ M, i.e. 
the level goes to the top of the gap, whereas in the limit 
of high R (large solitons) the behavior eo(M) ~ —M + a 
shows that the level goes to the bottom of the gap. The 
small M behavior, together with the reflexion property 
Eq. I[68p suggests a branch point behavior at the origin 







(79) 



which is no mystery, since it corresponds to a free par- 
ticle with the complex mass and with energy \/k 2 + ui 2 
at zero momentum. Unfortunately, in our study we only 
have access to the chiral soliton spectrum for real M, 
thus we do not possess the information on analyticity 
properties of the eigenvalues as functions of w, which is 
fully accessible only in exactly soluble problems. 

Indeed, due to the inherent numerical complications, 
finite-size discrete basis used, etc., such information 
would be difficult to obtain numerically in a reliable way. 
Thus, we proceed motivated by the real mass results and 
the analogy to a similar exactly-solvable model presented 
in Appendix [Dj In other words, our assumptions made 
in the general hedgehog case appear to be justified. 

Let us now review the calculation of the Dirac propa- 
gator for the standard case of a real mass. If we first com- 
pute the v integral in Eq. ([6"9"|) , we have to consider the 
7 contour in the complex-^ space which for real masses 
and hence real eigenvalues has the standard form of going 
below the real axis for states with energy below the gap 
and going above the real axis for states above the gap. 
This yields 



dv 



y(t-t') 



, 7 2m v - e n {M) 
[9(t - t')e(e n (M)) - 0{t> - t)9(-e n (M))] e -^(M)(t- t '), 



hence the positive (negative) energy states propagate for- 
ward (backward) in time. This is equivalent to a damped 
(exploding) imaginary-time behavior for t — t' — > ±«oo. 
This real mass result coincides with the standard one. If 
we change M the valence level may change sign (because 
it is equivalent to changing the soliton size R and the 
level dives into the sea according to (|78[)). Thus, accord- 
ing to Eq. ([50)) one keeps the contour. If we deformed the 
contour deciding that we occupy the negative eigenvalue 
then the imaginary-time behavior implying the standard 
particle (antiparticle) interpretation would be violated. 

Let us now analyze the case of a complex mass. For 
that purpose we go slightly off the real axis, u = M + iT, 
and expand 

e„(M + iT) = e n (M) + iTe' n {M) - \T 2 C(M) + . . .(81) 

The edges of the bound state gap and consequently the 
Dirac eigenvalues wander into the complex-^ plane mov- 
ing upwards or downwards, depending on the sign of 
e' n (M). The eigenvalues fulfilling e' n (M) = stay sta- 
tionary on the real axis. According to Fig. J2]), we have 
Cq(M) < for any M > 0. The question is whether or 
not we should also deform the contour 7 in the evaluation 
of the propagator ([6"9")l . On the one hand, the eigenvalues 
lying in the complex plane should not cross the contour, 
as this would lead to discontinuities, thus some deforma- 
tion must be done. On the other hand, if e' (M) = we 
move along the real axis in the positive direction by the 
amount — |r 2 £Q(M), Eq. ([8Tj) . hence crossing the con- 
tour only in the case where also eo(M) = 0. According 
to Fig. © this situation never happens unless M = 
or M is much larger than the position of the maximum 
(which, as we will see, is of no relevance to our construc- 
tion). Hence, we deform the contour in such a way that 
it continues to pass through the v — point, but never 
allows a complex eigenvalue to intersect. Such a contour 
yields 



2m' v — e n (ui) 



- 6{t' -t)0(-Ree„(w))]e-"" (w)( *- t ' ) . (82) 



We note that the large imaginary time evolution is 
damped, as we might have expected. In the positive 
t—t' > branch, where a baryon (and not an antibaryon) 
propagates, we are thus left with the integral 



duj P (Lo)9(Ree n (Lo))e- Te " {uj) . 



(83) 



for the imaginary time r — > 00. Note that we only have 
the positive section of the contour, as implied by the con- 
dition Ree„(w) > 0. Recall that we wish to occupy the 
lowest + orbit with the eigenvalue eo(w). We proceed by 
using the saddle-point method. First we have to locate 
stationary points of eo(w), i.e. the points uj m where 



e' (uj m ) = 0. 



(84) 
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FIG. 3: The deformed contour C' in the complex tj-plane for 
the saddle-point evaluation of the valence contribution. 



From Fig. [3 we find, that for any value of R there exists 
a saddle located at the maximum of the curve eo(M). 
We denote the position of the saddle as Mo 11 . We call 
Mq the saddle mass. As we move left from the saddle, 
decreasing M , we enter at some critical value M c into the 
upper continuum. On the basis of the toy model results 
of Appendix [Dl we expect that eo(w) has a branch cut 
at lo = M c running downwards. We observe from Fig. [2] 
that the saddles are located right of M c for any value of 
the soliton size R. 

In order to compute the integral in Eq. ([55)1 we deform 
the original contour C used in the vacuum sector into 
the contour C shown in Fig. [3) It contains the segment 
C — C (dashed line) , which runs across u> = Mq parallel 
to the imaginary axis. Note that it is not necessary to 
change the path globally but only in the vicinity of the 
stationary point. Making the change of variables lo — 
Mo + iT, where M is the saddle mass, e'(M ) = we get 

f dcjp(ij)9(Ree n (u;))e- Te " {uj) 
Jc-c 

/oo 
idr P (Mo)e- Te ^e^ M ^ r2 
-oo 

V - e oWo)r 



A priori we should admit multiple saddles uj m . We choose the 
branch, denoted u>o, with the lowest possible real part of the 
eigenvalue. In addition, the complex reflexion property implies 
that if luq is a saddle, then u)^ is also a saddle. As a consequence, 
for complex saddles the large Euclidean time behavior would 
not be purely exponential, but also some oscillations would set 
in. This obviously contradicts the spectral decomposition 160 1 , 
which allows only for real stationary points. Therefore, although 
complex solutions may a priori exist, they should be considered 
a spurious result of the model. Similar problems occur also in 
the local models fsstl . 



= Z e' T < M °'> (85) 

where the wave function renormalization factor Zq has 
been included. Finally, we compare Eq. ([55)) to Eq. ([tjQ]) 
and read off the valence contribution to the energy as 

E val = JV c e (M ). (86) 

The necessary conditions are e (M ) — 0, e (M ) > 0, 
and e '(Mo) < 0, all satisfied in our case, as shown in 
Fig. [2 In addition, the solution is only admissible with 
the conditions M c < Mq < My/2, which guarantee a 
real residue. As we have already discussed, M c < Mo 
holds for all values of R. For M > My/2 we would have 
a complex residue, contradicting the spectral decompo- 
sition of the propagator, Eq. ([60]) . According to Fig. [2[ 
for the ansatz ([T2"[) this occurs for R < 0.6 fm, where no 
bound-state nucleon can be constructed. Thus, in the 
large imaginary time limit t — > — ioo and t' — > +ioo we 
get ( see e.g. Eq. ([55]) ) 

S(x\ x) -> ZoMx)*o(x')e- i{t - t,)eo , (87) 

where 'i'o is the valence Dirac spinor. 

To summarize, our prescription for the valence orbit 
consists of the following very simple steps: 

1. Obtain (numerically) the spectrum for the valence 
orbital as a function of the real positive mass M. 

2. Look for the saddle mass, i.e., a maximum with 
respect to M, with the help of conditions e (Mo) = 
and e '(Mo) < 0. From analyticity, this is the 
saddle point which corresponds to the saddle mass 
M . 

3. If M c < Mo < My/2, the valence contribution to 
the soliton mass is N c eo(AIo), otherwise there is no 
baryon state in the model. 

There is an essential difference between the standard va- 
lence prescription based on the Fock-space decomposi- 
tion, used in local quark models, and our prescription for 
SQM presented above. In the first case the size of the 
profile, R, and the quark mass, M, can be fixed inde- 
pendently in the Dirac operator, whereas in the present 
case we actually find a correlation between them. In 
the quark propagator we have a superposition of masses 
to which are integrated over 12 , thus no additional inde- 
pendent scale is present in the problem. An immediate 
consequence can be seen in Fig. 21 where we compare the 
valence eigenvalue for the profile function ([5D[) along the 
manifold where e (M) = (for a given R) to the fixed- 
mass result of the local chiral quark model, obtained with 



This reminds of the familiar distinction between the phase and 
group velocities of wave packets made out of a continuous super- 
position of plane waves. 
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Mq = 300 MeV. From the scaled Dirac equation J76 
and our valence prescription it follows that 



e = k/R, 



(88) 



where k is a constant depending on the particular profile. 
In particular, for R — > in SQM one has eo — > 1/R (a 
signature of a "repulsive force"), as opposed to the be- 
havior eo - ¥ Mq of a free particle at rest in the case of 
the fixed-mass model. This means in practice that the 
bound state never becomes unbound in the presence of a 
chiral field background (see the discussion of Sect. IIII F|> 
and reflects the absence of on-shell quarks in the vacuum 
(analytic confinement) in SQM. Another feature of SQM 
with our valence prescription is that the valence eigen- 
value never dives into the negative part of the spectrum. 
The prescription for fixed- mass models allows the valence 
contribution to become negative, and this happens with 
a finite slope which originates from non-analyticity. This 
entering of the valence level into the negative part of the 
spectrum was interpreted within the chiral quark soli- 
ton model as entering the "Skyrme" model regime. In 
our approach such a regime never arises. However, the 
1/R dependence of the valence contribution behaves very 
much like the Skyrme stabilizing term. 
Similarly, the saddle mass scales as 



M = k'/R. 



(89) 



Numerically, for the exponential profile ([72]) one gets k = 
123 MeV fm and k' = 151 MeV fm, for the linear 
profile (73]) fc = 234 MeV fm and k' = 257 MeV fm, 
and for the arctan profile (f74"]) k = 100 MeV fm and 
k' = 119 MeV fm. 



D. The soliton baryon density and the baryon 
number 



A crucial point is to show that our soliton corresponds 
to a system with baryon number equal to one. This re- 
quires the separation between the valence and sea con- 
tributions to the baryon number. To analyze this point 
we have to compute the following time ordered product, 

U%(x,x',y) = (0\T{B(x)q(y)^q(y)B(x')}\0) 



(90) 



with B(x) denoting the interpolating baryonic operator, 
Eq. (|58[) and H v B {x : x') the baryon correlator (j5"7[) in the 
presence of a external vector field v fl (y) for which the 
result (|6T|) follows if the external source is included in 
the propagator 

SZa>M= I <M")<*l(i#-^-«tf 6 )~V>-(91) 
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FIG. 4: Solid line: the valence eigenvalue eo and the saddle 
mass, Mo, for the + state (maxima of Fig. [5J, plotted as a 
function of the soliton size R. Dotted line: the + eigenvalue 
e at the fixed value of M = 300 MeV, as used in traditional 
fixed-mass models. We note that the SQM prescription has 
eo = const /R, whereas e — > const at R — » and becomes 
negative at R — > oo. 



For the open line propagator the functional derivative 
can be readily evaluated yielding 



5S v (x',x) 



8v^{y) 



du>p(uj) 



x (x\ 



1 



-\y)r(y\- 



i 



*o(yh^o(y)S(x,x') 



(92) 



where in the last line the limit t — *■ —zoo and t' — ► +ioo 
has been taken along the line of reasoning developed in 
Sec. IIII CI (see Eq. |8"T|) . The determinant contribution to 
the baryon current can be deduced from Eq. (|5ip . Col- 
lecting all results we obtain for the three-point correlator 
in the asymptotic limit the factorized form 



n»(x,x',y) -> B»(y)n B (x,x>) 
with the total baryon current given by 



(93) 



B»{x) = * (x)j^ (x) 



diup(uj)tT 



ujU 5 ' 



(94) 



which is our result for the baryon current and in gen- 
eral for any observable based on a one-body operator. In 
Fig.[5]the sea and the valence contributions to the baryon 
density for several soliton radii in the SQM are displayed. 
As we see, the Dirac sea contribution gives a vanish- 
ing contribution to the baryon number; only the valence 
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FIG. 5: Sea and valence contribution to the radial baryon 
density for several soliton radii in the SQM (the valence parts 
have positive sign at all r). 



quarks contribute to the total normalization. This is con- 
sistent with our finding that the valence level never dives 
into the negative energy region. 

It is worth mentioning that the previous calculation 
can be extended mutatis mutandis to any one body ob- 
servable and in particular to the energy-momentum ten- 
sor. In the following section we analyze the sea contri- 
bution as it arises from Eq. (|60p onwards. A non triv- 
ial check which results from the one body consistency 
condition is that such a calculation coincides with the 
one based on the energy-momentum tensor in the soliton 
background. 



E. The sea contribution 



fined in Eq. ([65)1 and the frequency integral has been 
carried out. In the last line Eq. (|6"T)) has been used and 
summation i runs over negative energy (sea) states only. 
A full-fledged evaluation of this quantity is presented in 
Sect. lIVl To understand the general trend we will present 
an estimate based on the so-called two point function ap- 
proximation (TPA) proposed orig inally in Ref. [20l | and 
further exploited in Ref. ffdl , |76| . This approximation 
has the virtue of reproducing both the limit of large and 
small soliton sizes and is based in the identity for the 
normal parity contribution to the action 



dwp(uj)Tr\og (i^-tuU 5 



c 



= - I dup(Lj)Ti log (d 2 + iusfiU 5 + lu 2 ) 



(96) 



(deduced by commuting the 75 matrix across the Dirac 
operator) and further expanding the logarithm to second 
order in the field U 5 , yielding 



E. 



TPA 



-iN c 
AT 

-iN c 
T 



duip(ui)Tr ■ 



c 



Q2 



T iu@U'- 



dujp(uj) 



c 



d A q 

W) 4 



q 2 I(q 2 ^)iU 2 (\U{q)\ 2 ) 
(97) 



where the functional trace and Dirac traces have been 
evaluated, (.) indicates the isospin trace, and U(q) is the 
Fourier transform of the chiral field U (x) (see Appendix 
[5] for notation). The one- loop integral I(q 2 ,oj) is intro- 
duced and calculated in Appendix [B] A straightforward 
calculation in the meson-dominance version of the SQM 
for static fields yields 



E c 



where 



TPA 

sea 



1 



-f 



d 3 q 
(2tt) 3 



M 2 



M 2 ' 



(98) 



To identify the sea contribution we compare, as done 
for the valence case, the spectral decomposition of the 
two point correlator a large Euclidean times, Eq. (f6"0")) . 
with the path integral representation, Eq. (fBTj) . Using 
Eq. (|6"3")) for stationary configurations, the Dirac sea con- 
tribution to the energy is obtained by the spectral inte- 
gration along the original contour C of Fig. [TJ 



E H , 



-N c / dcjp{uj)Tr log (i$ - ujU 5 ) 
1 Jc 



-N c / dLop(u)\ y/e n (u 

J C 2 n 

N c / dujp(uj) ^ £i(w), 
Jc .TZ„ 



(95) 



where T is the time and n indicates all (positive- and 
negative-energy) states of the Dirac Hamiltonian H de- 



U(q) = / d A xU{x)e l 



(99) 



The factor appearing in Eq. (|98|) is the pion form fac- 
tor, Eq. (fl"8"|) in the space-like region, q 2 = — q 2 < 0. 
Actually, this is a general feature, static soliton profiles 
probe the Euclidean region of mesonic correlation func- 
tions. At large values of the soliton size R, small values 
of q dominate in Eq. (f9"8")l and one gets 



(100) 



E ~ a/ 2 i?, 



where a is a numerical constant. For the exponential pro- 
file ([72]) one gets numerically a = 30.99. In the opposite 
limit of small soliton sizes, large q values dominate and 
Eq. ([98} yields 



E ~ bf 2 R(M v R) 2 



(101) 



16 



with b = 28.11. This behavior is different from the 
R 3 log R short distance behavior documented in Ref. [2(| 
for the proper-time regularized fermion determinant; it is 
also free of the R log R behavior reported in Ref. [77], [78| 
for the renormalized sea energy, which generated a vac- 
uum Landau instability. 



F. Existence of absolute minima 

The total soliton energy is the valence plus the sea 
contribution as 
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(102) 
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At small radii the valence contribution (|88|) dominates 
and E ~ l/R, while at large radii the Dirac sea contri- 
bution dominates and E ~ R. Since the function E(R) is 
continuous, on these simple grounds we prove to have a 
minimum. The behavior is illustrated in Fig. ([6]), where 
we show the energy dependence as a function of the pro- 
file size R for the profile ([72"]) . For this restricted con- 
figuration one clearly sees the occurrence of an absolute 
minimum as a consequence of the fact that the valence 
quarks never become unbound. In contrast, in the stan- 
dard approach to chiral quark solitons the valence quarks 
become unbound at small R and one has instead a local 
minimum, which becomes unstable when the total soliton 
energy exceeds that of N c free quarks at rest, E > N C M. 
This instability has been an obvious cause of concern 
since this situation corresponds to bound but not con- 
fined solitons. In practice, it would not be a problem if 
valence states where deeply bound. However, most calcu- 
lations have produced solutions which lie on the unstable 
branch. The fact that our solutions correspond to an ab- 
solutely stable state is a remarkable property of SQM 
together with our construction of the valence state. Re- 
sults for other profiles of the chiral field are qualitatively 
similar to the case of Fig. [6] 

Another interesting feature that can be seen from a 
direct comparison of Figs. [6] and [4] is that the minimum 
takes place at a soliton size R where the valence state 
for our model and the fixed-mass models produce a simi- 
lar dependence. This suggests that the constituent fixed 
mass in soliton models may indeed correspond to a given 
saddle mass in the spectral construction and that the 
shallow (and unstable) minima found there may indeed 
be identified as the absolute and stable minima obtained 
here. 



IV. RESULTS FOR THE SELF-CONSISTENT 
SOLITON 

In this section we describe the numerical properties 
of self-consistent chiral solitons in SQM, as well as the 
corresponding observables. 



FIG. 6: Total energy of the soliton (solid line) and its valence 
and sea components for the profile (|72[) . plotted as functions 
of R. An absolute minimum exists. 



A. Evaluation of observables 

As is well known, observables in hedgehog solitons fall 
into two categories: independent of and dependent on 
cranking [7^ |. The second case is much more compli- 
cated, hence in this paper we deal only with the quanti- 
ties that do not involve cranking. Observables indepen- 
dent of cranking can be written in terms of single spectral 
sums of the form 

A = ^Ai, (103) 

i 

where A is the value of the observable in the soliton, and 
Ai is the single-particle contribution which includes a va- 
lence orbit and the Dirac sea. In the spectral approach 
this is generalized, in the sense that an extra integration 
over the spectral variable u> is present. With our method 
of treating the valenceness we also have a separated con- 
tribution from N c occupied valence states, 

A = A va \ + ^4 S ea, 
^val = N c A(u> va \), 

A sca = / dup{oj)A sca {u), (104) 
Jc 

where val and sea indicate valence and sea, w va i = Ma 
is the saddle mass, p(u>) is the spectral function con- 
sisting of the vector and scalar parts given in Eq. (|20[) . 
and A sca (uj) = Z) iesoa ^i(w) with A t (oj) denoting the in- 
dependent single-particle value of the observable in the 
orbit i. Due to the hedgehog symmetries for the ob- 
servables in question one may replace the sum over the 
negative energy states, X^igsca' w ^h 1/2 of the sum over 
all states, \ J2nea\i- 

The contour C is invariant under the transformation 
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see Fig. [T] Therefore 



A 

1 „ 



dujp(Lo)A sea (uj) - / du>p(-u>)A sea ,{-u)) 
uc Jc 

. A ae Ju>) + A ae J-Cj) 
dwpy (W ) 



du>ps(u> 



A sea (u>) - A sea (-uj) 



(105) 



where pv and are the odd and even parts of the spec- 
tral function p, respectively, see Eq. (|20|) . As we have 
stressed before, the contour C in Fig. [1] is complex. This 
is a complication, since in a numerical calculation we ob- 
viously do not have access to A(u>) in the complex ui 
plane. For that reason we use a method which allows 
us to carry on the calculation along the real axis in the 
u) plane, as explained in Appendix [F] According to the 
formalism of Sect. Mil the the energy of the soliton is 



E — N c E va i 
Eys(u) 

><E 



du 2 



[E v {u) + E s {u)] 



u=l/M? 



2 / cLudisc[py s(w)] 

Jl/(2y/u) 

6i(u>, m) ± e,;(— u>, m) 



(106) 



2 / dudisc[py s{u)} 

'l/(V5) 

ei(w,m) ± ej(o;, -m) 



E 



With the help of Eq. (|F15|) we may also write 
Ev,s{u) - 

In the chiral limit 
E s (u) = 0. 



(107) 



2 / _ disc[py(o;)]^e 4 (w,0), 



(108) 



All observables not involving cranking can be evaluated 
with the technique described above and with the help 
of "standard" formulas at each ui. The derivatives with 
respect to u are carried out numerically. 



B. Self-consistent equations 

Now we present some details of our self-consistent pro- 
cedure. We use the linear representation of the hedgehog 
field, U — s + ir-rp, and impose the nonlinear constraint 
s 2 + p 2 = 1 by renormalizing the fields after each nu- 
merical iteration. The Euler-Lagrange equations for the 
radial s and p fields are obtained from Eq. (|G1|) and have 
the form 

^ = N c Mx,M )Mx,M ) 
9 



c 



p(r) 



i 

N c Mx, M )i l5 T ■ rfa{x, M a ) (109) 
dupiui)^ i->i{x,Lo)i~j b T ■ fij}i(x,u)), 



c 



where g 2 is treated as a Lagrange multiplier. The first 
terms on the r.h.s. are the valence contributions, eval- 
uated according to the SQM prescription. The second 
terms are the Dirac sea contributions, where the sum 
over the negative-energy states i is carried out, or, equiv- 
alently, it can be replaced by 1/2 of the sum over all 
states.. With the help of formulas (|F10IF11|) the spectral 
integration in (|109[) can be performed as a real- valued in- 
tegral. The spinors tpi(x,uj) are obtained by solving the 
Dirac equation in the background of the fields s(r) and 
p{r) at a given value of ui. 

The code used to find numerically the self consistent 
solutions is a modification of the method used in solving 
different versions of chiral-quark models. The quark or- 
bits are calculated by diagonalizing H for each value of 
lo in the discrete Kahana-Ripka basis [l6j]. The Euler- 
Lagrange equations are solved by iteration. The numeri- 
cal effort involved in the calculation is similar to that in 
the case of solitons in non-local models (48l HH . 



C. The self-consistent solution 

For the sake of obtaining physical properties of the 
soliton, we use two versions of the model. The first 
one (model I) is just the standard meson-dominance 
SQM, with the single vector meson mass, My — m p — 
769 MeV. This model (in the chiral limit) has only 
one scale, thus all observables are proportional to the 
value of My in appropriate power. In particular, f 2 — 
M v /(8ir 2 ). The second model (model II) includes also 
the excited p state, p'(1465). The spectral function is 
taken to be the weighted sum of Eq. ([2H)l . containing 
90% of the ground-state p, and 10% of excited p. All 
quantities are distributive over the spectral density, e.g. 
f 2 = 0.9m 2 p /{8n 2 ) + 0.1m' 2 /(8Tr 2 ). Model II contains two 
scales, and produces somewhat heavier and more com- 
pact solitons, as expected on simple scaling grounds. We 
work with the physical pion, — 139.6 MeV, and the 
current quark mass is adjusted such that the GMOR re- 
lation is fulfilled, m(qq) = m^f 2 . 

The results of the self-consistent calculation for the 
case with physical pion mass are displayed in Table U 
and in Fig. [7] The chiral corrections are small and do 
not alter our basic conclusions. We note that for both 
models the value of the saddle mass Mo is well below the 
critical value \My discussed in Sect. Mil 

Comparing the soliton energy to the experimental nu- 
cleon mass we should be aware that our soliton is a 
mean-field solution with grand spin and should not 
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r[fm] 

FIG. 7: The s = cos 8 and p = sin 8 fields as a function of the 
radius r for the two self-consistent solutions of Table [T] the 
solid lines for the model I and the dashed lines for model II. 



be identified with the nucleon, but rather with the av- 
erage of the mass of the nucleon and the A(1232) iso- 
bar \{M N + Ma) = 1174 MeV. Moreover, quantization 
of the collective coordinates, or projection of the soliton 
wave function on the subspace with good quantum num- 
bers [13, HH , reduces its energy by eliminating the spuri- 
ous rotational and translational energy. We note that the 
contribution to the energy arising from the scalar spec- 
tral density, E scg ^s, is much smaller than the vector part, 
E>sea.,v- Furthermore, we notice that the model gives ap- 
proximate equipartition of energy between each valence 
quark and the Dirac sea, in agreement to the 1 / R behav- 
ior of the valence contribution and the approximate ~ R 
dependence of the sea part. 



model I model II 


U [MeV] 


86.5 


97.3 


m,r [MeV] 


139.6 


139.6 


(~(qq)) 1/3 [MeV] 


243 


243 


m [MeV] 


5.04 


6.37 


M [MeV] 


267 


304 


£ [MeV] 


233 


263 


£ sea ,v [MeV] 


285 


351 


£ S ea,s [MeV] 


36 


35 


Stotai [MeV] 


1019 


1174 



TABLE I: Model parameters and the soliton energy in the 
self-consistent calculation with the physical pion mass. 



D. Results for observables 



model I model II experiment 



\/{7- 2 )i=0,val [H 

y/ <r 2 )i=o,totai [fm] 


1.23 
1.24 


1.07 
1.08 


0.79 


ffA(val) 
gA (total) 


0.79 
0.93 


0.79 
0.95 


1.26 


^/=i(val) 
(total) 


2.40 
2.96 


2.07 
2.67 


4.71 


ct^n [MeV] 


36 


33 


45 ±8 



TABLE II: Properties of the self-consistent solution of Table|T] 
the rms isoscalar charge radius, the axial- vector charge gA, the 
isovector magnetic moment (j,i=i, and the sigma commutator 
ovn- The valence contribution (val) is given separately. The 
contribution from the sea (not given explicitly) is dominated 
by the vector part, similarly as in the case of the energy. The 
experimental value of ct^n is taken from Ref. [s3 ] ■ 

cranking. Inclusion of cranking, which requires a linear- 
response calculation, is numerically involved and is out- 
side of the scope of this study. In Table[TT]we display some 
characteristic observables of the self consistent solution 
for the choice of parameters of of Table Q] Firstly, we 
note that the isoscalar rms charge radius of the soliton is 
too large as compared to the experiment. This quantity 
is dominated by the valence contribution, as can also be 
seen from Fig. [8] The quark spinors of the valence wave 
function exhibit a long tail which is a feature related to 
our prescription for constructing the valence orbit. The 
large r behavior is of the form exp(— \/Mq — e§ r), and 
in our case the value of k = \j M§ — is small, with 
k = 1/(1.5 fm) in model I and n — 1/ (1.3 fm) in model II. 
The large soliton radius can perhaps be reduced by sub- 
tracting the spurious center-of-mass motion of the soli- 
ton. Including higher vector mesons in the model tends to 
reduce the soliton size, as expected by scaling arguments, 
however, on the expense of increasing f v and -Etotai- 

In the evaluation 13 of gA and the isovector magnetic 
moment we use the semiclassical projection coeffi- 

cients as obtained in the large- N c limit [79[ . Both values 
are lower than the corresponding experimental values. 
But since our solution is dominated by the valence state 
we could as well, instead of the semiclassical projection, 
use the quark model wave functions for evaluation of the 
valence parts of observables. Then the expectation value 
of the or operator for the proton and the neutron yields 
the coefficient | instead of | . That would result in higher 
values, in better agreement with the experiment. The 
value of the sigma commutator, which in our model is 
equal to 

t.n = -mlfl f d 3 r(s(r) - 1) (110) 



As discussed in Sec. IIV1 in the present approach we 
deal only with those observables that do not involve 



Explicit expressions for these quantities may be found, e.g., in 
Ref. [H,|27|. 
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FIG. 8: The valence and the sea contributions to the baryon 
density (multiplied by 4-kt 2 ) for the self-consistent solu- 
tions:the solid lines for model I and the dashed lines for model 
II. 

is reasonably well reproduced, assuming values somewhat 
smaller than the experimental number, showing that the 
spatial extent of our chiral profile is not too large. 

The results for the observables turn out to be quite 
similar to those obtained in other chiral models. In par- 
ticular, they are strikingly close to the predictions of the 
model with the non-local regulators [48|, Eg] • 



V. CONCLUSION 

For many years it has been accepted that baryons arise 
as solitons of a chiral Lagrangean in the large- iV c -limit 
of QCD. However, the precise realization of this attrac- 
tive and fruitful idea is still unclear. Many model La- 
grangeans have been proposed, emphasizing different as- 
pects of the problem and incorporating as many known 
features of the underlying quark-gluon dynamics as pos- 
sible. While on the one hand the Skyrme soliton models 
incorporate confinement but not the spin-1/2 partons, 
on the other hand, chiral quark soliton models account 
for spin-1/2 partons but lack the confinement property. 
This unsatisfactory situation has been calling for im- 
provement. As a first and hopefully useful step, we have 
considered a model where quark poles on the real axis are 
absent, be it in the vacuum or in the soliton. Unfortu- 
nately, confinement does not hold since basic large- N c re- 
quirements on the analytic structure of meson correlators 
are violated, similarly to other effective quark models. 
On the other hand the Spectral Quark Model provides 
a framework with spin-1/2 partons in the vacuum which 
simultaneously allows for baryonic soliton solutions cor- 



responding to absolute minima of the action and hence 
cannot decay into three free quarks, unlike previous local 
chiral quark soliton models. This is due to the rather 
peculiar features of the model characterized by a contin- 
uous superposition of masses in the complex plane with 
a suitable spectral function. As a necessary complemen- 
tary study we have examined an instance of an involved 
complex-mass relativistic system. 

The solitons we find fit nicely into the phcnomcnolog- 
ical expectations of more standard chiral quark solitons 
with full inclusion of the polarized Dirac sea. In contrast 
to constituent chiral quark models where one is allowed to 
tune the constituent quark mass as a free parameter, the 
Spectral Quark Model does not have this freedom. In our 
case the vector meson mass determines one unique solu- 
tion in the chiral limit and thus all properties scale with 
this mass. Actually, due to the superposition of com- 
plex spectral masses u>, in the soliton calculation there 
appears a saddle mass, Mq, defined as a stationary point 
of the valence quark eigenvalue as a function of lu in a 
fixed soliton profile. In the chiral limit Mq also scales 
with the vector meson mass and thus is a fixed num- 
ber. To some extent the saddle mass behaves as a con- 
stituent quark mass determined uniquely by the soliton, 
and its numerical value comes out within the expected 
range Mq ~ 300MeV for the typical solitons minimizing 
the total energy. Moreover, the close resemblance of the 
valence contribution to the energy around the minimum 
suggests that the phenomenology and results of previ- 
ous calculations within constituent chiral quark models 
with a constant mass are reproduced, at least for the 
low-lying baryon states. The new aspect unveiled by our 
calculation is the possibility of studying the monopole vi- 
brations and their quantization which are the traditional 
candidate for describing excited baryonic states such as, 
e.g., the Roper resonance. While studies of this sort are 
routinely carried out in the Skyrme model, the lack of 
confinement in the traditional chiral quark soliton model 
has prevented, as a matter of principle, calculations of 
excited states. 

One distinguishing feature of the Spectral Quark 
Model construction is the verification of an important 
consistency relation involving quark two-point functions 
which were not fulfilled for fixed constituent mass chiral 
quark models. This feature, in addition to the uniform 
treatment of regularization removes theoretical doubts 
on the proper computation of the high-energy processes 
as regards their partonic interpretation and normaliza- 
tion. As we have mentioned, one motivation for using 
chiral quark soliton models is the spin 1/2 nature of the 
constituents, which makes the calculation of partonic dis- 
tributions in the nucleon possible as a matter of principle. 
The results found in the present paper pave the way for 
such a calculation in the soliton picture where the in- 
terplay between the analyticity enforced by the lack of 
on-shell quarks and the chiral symmetry may simultane- 
ously be tested for the low-lying baryon states. 
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APPENDIX A: THE EFFECTIVE ACTION TO 
SECOND ORDER IN THE FIELDS 

When using the path integral (JS]) with the effective ac- 
tion (fl7J| it is enough to expand the action to second order 
in the external sources and the dynamical pion field, 



U 



1 + jT- 7?7 5 



1 



(Al) 



This can be done by noting that the Dirac operator in 
Eq. (fTTj) can be decomposed as a free part plus a pertur- 
bation, 



iD = i$ - oj - V, 



(A2) 



where 



V = UJ \^jf ■ 7T75 - 7? • 7T + . 

+ m - {i> + 4l5 - s - ij 5 p) 



(A3) 



Note that both the free propagator and the potential V 
may depend on the spectral mass uj. Then the perturba- 
tive expansion of the fermion determinant can be readily 
done, yielding 



iT = N c 



dujp(uj)Ti log(z^) — uj) 



N 



f (kj P {u>) jr -ti 

Jc ±i« 



-v 



(A4) 



which can be classified according to the number of exter- 
nal sources as well as dynamical fields. 

The zeroth order contribution yields the vacuum en- 
ergy, 



ir (o) = N c I dujp(w)Trlog(i@-w) 
Jc 



N c N f 
= I d 4 x i 



d^x I dujp{uj) 
c 



d 4 p 



trlog(/5 — uS) 



N e N f 



d 4 p 



trlog(^)L (A5) 



where the vacuum energy density, e, defined relative to 
the case of massless quarks, has been explicitly separated. 
This quantity was computed in Ref. [53| and in Sec. Ill CI 
of the present paper through the use of the energy mo- 
mentum tensor. It is instructive to check the calculation 
in the effective action formalism. We use the reflection 
symmetry p — > — p and the standard identity for matrices, 
trlogA+trlogB = logdet A+logdetB = logdet(AB) = 
trlog(/l£>) for the commuting Dirac matrices A = j> — uj 
and B = p 1 + lo (in four dimensions det A = det(— A)). 
The vacuum energy density can be written as one half of 
any of these contributions. The Dirac trace becomes triv- 
ial and after interchanging the order in the momentum 
and spectral integrals one gets 



-2iN r N 



d 4 p 



(2tt)4 

duip{uj) [log(p 2 - uj 2 ) - \og{p 2 )} 

To evaluate first the w-integral we write 

1 dk 2 



\og(p' - OJ J = 
Thus, we readily get from direct application of Eq 
d 4 P 



(A6) 



(A7) 



-2iN c N 



iN r N 



(2tt)4 
d 4 p 



dk 1 



p 2 A( P 2 



A(k 2 ) 
1]. 



1 



(A8) 



after going to the Euclidean space and integrating by 
parts. The surface term has been discarded, as the condi- 
tion linis^+oo s 3 A(—s) = 0, holds in the meson dominance 
realization Eq. (|2"Tj) . Thus the final result is finite, 



N c NfMy 

1927T 2 



-(0.2GeV) 4 (N f =3) 



(A9) 



in agreement with Eq. (f2"7) based on the energy- 
momentum tensor. 

Returning to Eq. (|A4j) . the first non- vanishing correc- 
tion involves only the scalar source, yielding 



Nr. 



4N C 



diop(uj)Tr 
d 4 p 



1 



(2n) 4 



duo- 



1 - UJ 
p(u>) 



(A10) 



c 



P 



d x(s(x)) 



((.) means the flavor trace), whence the quark condensate 
may be obtained as 



Nf(qq) = —AiN c Nf 



d 4 p 



B(P 2 ), 



(All) 



which becomes an identity according to Eq. (|2"Tj) . 
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Finally, the second order contribution can be written 



as 



x Tr 



„ 1 i 
Ta~, ; T B - 



p 1 — $ — uj j) — uj 



.(A14) 



r (2) = T ss + r PP + T vv + T AA + r Wff 

+ r nP + r vA + r PA . (A12) 



The Fourier-transformed fields are defined through the 
The calculation is straightforward. For the bilincars in relation 
the fields, generically written as (f(x) = ip A (x)T A with A 
a denoting Lorentz-flavor index, we have 



1 



N c / dujp(uj)Ti 



2 

; f d 4 q 



up — UJ 

2\ 



2 J (2tt)' 



■tpA(q)VB(-q)K AB (q ), 



(A13) 



where 



(-i)K AB (q) = -N C J dup{u,) J 



(p(x) 



d 4 q 



(A15) 



The flavor trace is trivial, (A a A/,) = 8 a b/2. After using the 
vanishing condition of the positive moments, Eq. (|23[) . 
the calculation of the Dirac trace and momentum inte- 
gration yields 



(-i)K a ss(<l) = -Nc / dup(u) 



-*)K a A T(q) 



d 4 P 



Tr 



A a i X b i 



2 p-j-uj 2 p'-uj 



5 ab N c / dujp(uj) [2I(0,uj)uj 2 + I(q 2 ,uj)(4:Uj 2 - q 2 )] 
Jc 



(-i)K P b p(q) = -N c / dujp{uj) 



(-i)K^(q) 



d 4 p 



Tr 



' ■ ■•, *75 c, 

■ uj 2 p — uj 2 



5 ab N c / dujp{uj) [2/(0, uj)uj 2 - I{q 2 ,uj)q 2 
Jc 



N c / dujp(uj) 



\s ab N c i ,/ 



d 4 P 

W) 4 

q^q v 



Tr 



A a 



A" 



r J Jc 



i-i-w ^ 2 i-uj 2 

«2 



dujp{uj) 



= —N c / dujp(uj) 



d 4 p 

W) 4 



Tr 



41(0, 



A" 



24tt 2 



2I(q J ,uj)(2uj 2 + q / ) 



p— $ — uj 2 p — uj 2 



dujp{uj) 



<5 ab 4iV ( 



q^q v 



Q / Jc 
dujp{uj)I(q 2 ,uj)uj 2 , 



4/(0, uj)uj 2 - JL- + 2I(q 2 , uj){Auj 2 - q 2 ) 



H)K a A f(q) = -N c / dujp(uj) 



d 4 p 

W) 4 



Tr 



A a 
■7m75tt 



_. A^ 

f-i-uj 2 |5- w 175 2 



S^iNcq 11 / dujp(uj)ujl(q 2 ,uj), 
Jc 



(-i)K£(q) = ~N C dujp(cj) 



d 4 p 

W) 4 



Tr 



75 A : -7 5 A 6 



/ 2 - u; i>~i~u> f t-ujf 



-S ab 4N c % I dujp{uj)uj 2 I(q 2 ,uj), 



H)KTf(l) = 'Nc I dujp{uj) 



d 4 p 

W) 4 



Tr 



i A a i . h uj 

77 — 1 7^75^- 77 H5A -r 

p-£-w 2 p-uj f 



-6 ab 4iN c ±- / dujp(uj)I(q 2 ,uj)LJ 2 , 
J Jc 



(A16) 



(A17) 



(A18) 



(A19) 



(A20) 



(A21) 



(A22) 
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(-*)#&(?) = -# c / dwp{iu) 



d 4 p 



Tr 



6 ab 2N c - / dw/9(w) [2/(0, 
J Jc 



i X a i b ui 

I(q 2 ,uj)q 2 uj] . 



(A23) 



where the basic one-loop two-point integral I(q 2 ,uj) is 
defined in Eq. (|Bl|) . 

To compute the correlation functions as functional 
derivatives we must take into account the contributions 
from the pion pole. This is accomplished in a standard 
way by eliminating the pion field at the mean field level 
through the equations of motion, ST/5ir a (x) = 0, yield- 
ing at lowest order in the field 

v a {q)Kt( q ) + p a (q)K: b p (q) + a^(q)K^ b (q) = 

(A24) 

and reinserting the pion field into the effective ac- 
tion f^ 2 ). This contribution exactly cancels the non- 
transverse piece of the AA correlator, reproducing the 
gauge technique result of Eq. (|36|) and an additional pion 
pole contribution to the PP correlator yielding Eq. . 



APPENDIX B: USEFUL INTEGRALS 



The basic two-point integral is given by 1 



d 4 k 



1 



1 



1 



i J (2tt) 4 k 2 -uj 2 (q- k) 2 - uj 2 



1 



d A k 



dx 



i J (2-7r) 4 J Q [k 2 - uj 2 + q 2 x(l - x 
dx log [uj 2 - x(l - x)q 2 ] 



1 



(47T) 2 J Q 

For the meson dominance case one can use the moments 
method, based on expansion in powers of q 2 and resum- 
mation of the series. This method provides a useful check 
to our computations. However, for the purpose of illus- 
trating the issues of quark unitarity with our unconven- 
tional propagator it is enlightening to provide general 
formulas in terms of the vector and scalar components of 
the quark propagator. The u) integrals can be evaluated 
first by using Eqs. (f2Tj) . We get the useful identities, 



du> — 

c IP 



A'(p 2 ) = 

B'( P 2 ) = - I du 
Jc 

A{p 2 )+p 2 A'{ P 2 ) = 



p{uj) 



u>p(u>) 
[p 2 -w 2 ] 2 ' 

L0 2 p{<jj) 



duj 
c IP 
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There is a typo in Eq. (A4) of Ref. 



B(p 2 )+p 2 B'(p 2 ) = - I du 
Jc 



uj 3 p(lo) 



, (B2) 
c [p 2 - oj 2 \ 2 



The normalization of the pion field yields 

p-"v(£,4,[*«f)]-yg- 



(2tt) 4 dk 2 



(B3) 



We also have 



dujp(uj)I(q 2 , uj) 



ii dx . 

Jo J (2*r) 



d 4 k ., 



-A' {k 2 + x(l - x)q 2 ) , 



dujp{uj)I(q 2 ,co)ui 



c 



i I dx 



d^k 



B' (k 2 + x(l - x)q 2 ) , (B4) 



and so on. The integrals can be evaluated by pass- 
ing to the Euclidean space k 2 — > —k% and (— i)d 4 k — > 
-n 2 k 2 E dk 2 E — ir 2 sds. We may then shift the integration 
variable in order to get, e.g., 



(Bl) J(<1 2 



i dx 



d A k 

WT 4 



16tt 2 



F (k 2 + x(\ - x)q 2 ) (B5) 



sds I dxF (— s + x(l — x)q 



i r°° r 1 

= 16^Jo dSF( ~ S) Jo dx [S + ^-^W} + , 

where we have introduced the distribution [x\ + = x6{x). 
Let us assume that for definiteness q 2 < 0. Then, the 
argument of the step function vanishes for x < X- and 
x > x+, where 




(B6) 



We have 



f dx [S + x{l - x)q 2 ] = S 
Jo 



Thus 



{-q 2 -AS)Jl + ^{q 2 + AS). (B7) 



-iiUf (B8) 

r- ? 2 /4 / 7^ 

J dSF(-S)Jl + ^(q 2 + 4S)j. 
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This formula can be analyticly continued to any complex 
q, and the second integral becomes a line integral in the 
complex S'-plane. Clearly, if the function F(—S) has a 
cut at say S = —My/4, the result of the line integral 
becomes path dependent and will develop an imaginary 
part discontinuity for q 2 > My. This is the way how a 
quark propagator with no poles generates the unitarity 
cuts in the two point correlators. 
We list the final results: 
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APPENDIX C: QUANTUM-MECHANICAL 
EXAMPLES OF COMPLEX-MASS SYSTEMS 



Establishing analytic properties of the Dirac operator 
eigenvalues for a complex mass is an involved mathemat- 
ical problem. To gain some insight and develop some 
intuition, in this Appendix we consider a few cases of 
non-relativistic quantum mechanical system with com- 
plex coupling potentials. Our aim is to define what we 
mean by the energy eigenvalues for a complex coupling, 
assuming that we know the definition for the real cou- 
pling, and to determine its analytic properties in the 
complex plane. The most natural and obvious way to 
do so is in terms of analytic continuation in the coupling 
from the real case. 

As a first example let us consider the harmonic oscilla- 
tor which has the reduced potential U — 2mV — m 2 L0 2 x 2 
and the ground state energy is given by Eq(lo) — lo/2 for 
real lo > 0. The problem is obviously invariant under 
the change lo — > —lo, and we should write Eq(lo) — \lo\/2 
for real lo. Clearly, we have a branch cut at lo = as a 
function of the coupling lo 2 , hence the right way to write 
the energy for a complex coupling is Eq(lo) — Vlo 2 /2. 
Written this way the bound state wave function is given 
by tpa = Ce - ^ 1 ^ 2 . The analytic continuation of the 
decreasing exponential becomes an increasing exponen- 
tial on the second Riemann sheet, lo 2 — > e 2m Lo 2 . For 
negative lo 2 we have a negative energy. So, the energies 
on the first and second Riemann sheets differ only in the 
sign. 

As a second example let us consider the hydrogen atom 
where we have V — —Z/r with Z > and the energy 



is Eq(Z) = —Z 2 /2. This suggests an analytic behavior 
in Z, and in particular having a continued bound state 
solution for repulsive potentials [Z < 0). Again, the 
bound state function u{r) = re" Zr transforms into a 
positive exponential for negative Z . So, we have a cut at 
Z = 0. For negative Z we have a real energy but it is 
on the second Riemann sheet. In this case the energy on 
the first and the second Riemann sheets coincide. 

As a final example, more directly related to the more 
complicated case of the toy model for the Dirac equation 
in Appendix [D] , we analyze the complex square well po- 
tential for which we have U(r) = — U$Q(R — r), with Uq 
complex. For real Uo the s-wave bound state solution is 
given by (we take 2m = 1) 



u(r) =Q(R- r)Ae~ Kr + 8(r - R)B sin Kr (CI) 



where 



K=^U Q -E 



(C2) 



The continuity condition for the logarithmic derivative 
yields the bound state relation 



A' cot KR 



(C3) 



which defines an implicit function Eq(Uo) which we want 
to extend to complex Uo- For real Uo there is a critical 
value Uo, c = ( 7! '/2R) 2 above which the equation has real 
solutions. Actually, close to the threshold we have for 
Uo > Uo.c the ground state 



Eo(U ) 



1 



R 2 (Uo - U ,c) 



1 - -1 ) (CTb _ U , c f 



(C4) 



which suggests an analytic behavior of the energy. Again, 
it is the wave function at large distances where we see 
that k — > — k if we loop once about the critical point 
Uo.c = {k/2R) 2 in the complex Uq plane. In this partic- 
ular case, this corresponds to the well known fact that 
a square well potential with a subcritical coupling has 
a virtual state and no bound state. As Uo — > + one 
gets k — » — oo or Eq — > — oo on the second Riemann 
sheet. Similar features should appear also for excited 
states, E n , with the corresponding critical values. Thus, 
we may define the continuous function Eo(Uo) for all real 
values of Uq. For Uq > Uo.c it corresponds to a bound 
state, whereas for Uo < Uq iC it describes a virtual state. 

If we go now to the complex plane in Uq we get the 
implicit function defined through 



A{k,U) = 



(C5) 



If we have a solution A(ko,Uo) = and go close to it, 
then by Taylor expanding we get 

= A(k + Ak,U + AU) 

= A( k0 ,U ) + dAi b Uc) Ak + dA ^AU 



dk 



dU 



(C6) 
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This way we can define a one to one relation unless either 
derivative vanishes at a critical point 



dA(k c , U c ) 
dk 



= 



dA{k c ,U c ) 
dU 



= (C7) 



If this is the case we have a square root branch point 
assuming 



d 2 A(k c , U c ) 
dk 2 



In our case we have the critical points located at 



k = Kcot(KR) 
kR 



(C8) 



(C9) 
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— cot{KR) (CIO) 
K 



Combining both equations we get 
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R 



1 

E 2 



1 = KRcot(KR) U c 



1+x 2 



(Cll) 
(C12) 



with x n = 0, ±4.49, ±7.72, ±10.9, . . .. Moreover, we have 



d 2 A(k c ,U c ) _ R 3 U C 



dk 2 



R 2 U C 



1 



(C13) 



which diverges for U c — 1/R 2 . Note that the critical 
points are located in the second Ricmann sheet. They 
generate branch points of second order, i.e. looping twice 
around the point the function returns to its original value. 
The Riemann surface is obtained by joining all critical 
points with a line. Since they are infinitely many, the 
cut divides the complex U plane into two disjoint pieces. 
The complete and extensive study of Riemann sheets of 
this particular problem can be looked up at Ref. [83| . 

The main outcome is that the analytic structure of a 
the analyticly continued eigenvalues of a complex cou- 
pling potential can be determined by the study of the 
critical points which generate branch cuts. Otherwise, 
the function is analytic in the complex coupling plane. 
This situation appears also in Appendix [D] for the Dirac 
operator with a complex mass. 



APPENDIX D: ANALYTIC PROPERTIES OF 
THE VALENCE EIGENSTATE IN A TOY MODEL 

The discussion in Sect. IIIII made explicit use of the 
fact that an integration path in the complex mass plane 
of the Dirac Hamiltonian can be deformed without pinch- 
ing any singularities. Although we cannot prove this in 
general, the complex-mass coupling Dirac systems may 
have unusual properties. In this appendix we investigate 
a model baring similarity to the full chiral model where 
our assumptions are verified. The main issue is both 
to define the meaning of an eigenvalue as a function of 



the complex mass to as well as to determine its analytic 
properties in the complex-mass plane. To our knowl- 
edge this topic is not discussed in the literature at all. 
However, strong similarities are found with the analytic 
properties of the energy eigenvalues of complex poten- 
tials, which have been motivated in the context of opti- 
cal models [U, [85[ and for review purposes Appendix [C] 
includes some warm-up quantum-mechanical problems 
which may be helpful in the understanding of analyticity 
properties that arise in studies of our type. The pre- 
scription for valenceness derived in Sect. HTT1 holds under 
specific conditions. In the chiral soliton model we have 
no precise knowledge of the analytic properties, hence it 
is not mathematically proven that the prescription can 
actually be used. The model of this section is much sim- 
pler and solvable semi-analyticly. It shows that the de- 
sired analytic properties are fulfilled, providing support 
for the method in application to the chiral soliton model 
described in this paper. 

Let us consider the Dirac equation for the state with 
grand-spin and positive parity, G p = + , with the upper 
component u and the lower component v. It has the form 

v! = -uM sin 6 - v(e + M cos 6) 

v' = u(e-Mcos9)+v(-- + Msm6) (Dl) 

r 

with the usual boundary conditions for a normalizable 

state 

u'(Q) = 0, v{0) = 0, tt(oo) = 0, v(oo) = 0. (D2) 

Following Ref. we look for a solution with the profile 
function 



(r) =wQ(R-r), 



(D3) 



where Q(x) the standard Heaviside step function. The 
profile (|D3|) has the winding number (6(0)— 6(oo))/w = 1. 
The analytic solution is 

/„ , k 3 sinh nr 

u{r) = Nle(R-r)—- 

\ e + M r 

+ 0(r — R) 



M-e r 



v(r) = N (o(R. — r) — (nr cosYiinr) — sinh(Kx)) 



+ 6(r - R)K(Kr + 1) 



(D4) 



where K = \J M 2 — e 2 and N is the normalization fac- 
tor chosen such that A^ 2 J d 3 x(u 2 + v 2 ) = 1. Using the 
matching condition expressing the continuity of u(r) /v(r) 
at r = R we get the eigenvalue equation 

kR(M - e) + 2M 



KRcoth(nR) = 



(D5) 



e + M 

which can be written it terms of dimensionless variables 
x = kR and xq = MR as 



x coth(ar) 



x(xo ± yjx\ — x 2 ) + 2x 
±V x o -x 2 +x 



(D6) 
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e (R) /M 




where ± corresponds to the positive (negative) energy 
states, e > (e < ). The value x = is formally always 
a solution of Eq. (|D6|) , but it does not correspond to a 
zero momentum state eo = M since the solution is 



1 2 M 

ti(r) = —Q(R- r ) + 6(r - R) 

R 1 r 



(D7) 



v(r) = Q(R-r) 



2Mr 
~3W 



1 



Q(r-R)^, (D8) 



which are regular both at r — and r = oo but fail 
to fulfill the matching condition since v/u — 2Mr/2> for 
r > R whereas v/u — 2Mr for r > R. 

The numerical dependence of the eigenvalue as a func- 
tion of the soliton size for the + state in the toy model 
is depicted in Fig. [9] For large values of the soliton size 
at fixed mass we get 



eo(R) - ^ 



(D9) 



which means according to the scaling property that this 
behavior also holds true for large mass and fixed radius. 
The solution is always in the positive energy region and 
never dives into the sea (despite having topological num- 
ber 1 - this was actually the point of Ref. [86]). For 
R — > oo one gets e (R) — > + . The reason is related to 
the fact that always sm9 = for the profile, and we have 
a scalar coupling for which i?7o7s = and the 

spectrum is symmetric. 

The threshold value for having a bound state is 
(MR) 2 > 3/4. Note that since kR = x changes sign 
across the threshold value, the bound state becomes vir- 
tual (exponentially growing). Thus, the analytic contin- 
uation of the unbound valence state for R 2 M 2 < 3/4 is 
an exponentially growing state. 

Equation (|D6[) defines implicitly a relation between x 
and xq. Actually, we can compute the inverse function 



x A (coth(x) + 1)' 
4(x + 1) (xcoth(a;) - 



(D10) 



FIG. 10: Solid straight lines are e = ±M, the solid curve is 
the bound state, the dashed curve is the virtual state. Blobs 
indicate the position of the saddle, the emergence of a bound 
state, and the bifurcation point of two real virtual states into 
two complex states. For the bound state eoR approaches unity 
at MR — > 00, while for the virtual branch it approaches the 
negative continuum. 



For x > the function on the r.h.s. is a monotonously in- 
creasing function and thus the inverse is uniquely defined. 
The minimum takes place at x — and assumes the value 
of 3/4. Hence there are no positive energy bound state 
solutions for Xq < 3/4. For x^ = 3/4 we have a negative 
energy solution kR — x ~ —0.848 and e ~ — 0.175/i?, of 
no concern here. For 0.512 < x\ < 3/4 we have a positive 
energy state with k < 0. For 0.512 > x 2 , the value of k 
becomes complex. For Xq —> we get x — > — 00. 
In summary, for the real-mass case we have 

• For 0.866 < xq we have one bound state solution. 

• For 0.797 < xq < 0.866 we have one virtual state 
solution. 

• For 0.712 < xq < 0.797 we have two virtual states. 
For xq = 0.7968 one at x = -0.797 and the other 
one at x = —0.196, the higher one being a contin- 
uation of the previous case. 

• For xq < 0.716 the two solutions collide at x = 
—0.574 and bifurcate into the complex plane. 

So clearly for xq < 0.797 we have a cut in the k plane. 
The situation for the real mass case has been displayed 
in Fig. 1101 This plot should be compared with the expo- 
nential soliton profile of Fig. [2l 

We can now proceed to the complex plane in uj. The 
obvious way to define the eigenvalue eo(w) as a func- 
tion of the complex variable uj is by analytic continua- 
tion, since we want to follow the evolution of the state 
in the complex-mass plane. For illustration purposes we 
plot in Fig. [TT] the corresponding analyticly continued 
eigenvalues when uj = e la for several soliton sizes R. 



2G 



Im e / M 
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FIG. 11: analyticly continued complex valence eigenvalue 
eo(R,w) in the toy model plotted as a function of the com- 
plex mass u = e IQ for a in the first quadrant, < a < n/2. 
The axes represent Re(eo) and Im(en). Subsequent curves, 
from right to left, correspond to increasing values of R = 
2, 2.5, 3, 3.5, 3.9. The points on the real axis are for a — 0, 
and the value of a parameterizes the curves, which evolve 
away from the real axis. The other quadrants are obtained 
by reflection about the real axis. 



Equivalently, we consider the function of the variable Xq. 
Thus, we take the solutions which were defined for real 
M fulfilling the boundary conditions and arrive at the 
transcendental equation (|D6[) . which now defines implic- 
itly eo(w). 15 Thus, if we differentiate with respect to x 
the implicit function we get the equation for the critical 
points 



2t 4 + t* 



4i 2 



2t + 1 + (2t + 3)V1^<2 



0, (Dll) 



where for simplicity the variable t — x/xq has been in- 
troduced. This is an algebraic equation with solutions 
t = t n . We find numerically 

t n =x/x = -1.2741 iiO.1174, 0.9092, 

1.3898 ±i0.4832,..., (D12) 

which forms a bundle of solutions 



t n x coth(i„x ) 



X t n (l-y/l-tl)+2t r . 

v / i^+ 1 



(D13) 



We see that for any value of t n there are infinitely many 
solutions for xq. Roughly, they are located along a 
straight line. For instance, for x/xq — 0.9092 the critical 
points in xq are located at about xq — 0.5 + imr for large- 
TV. So, the whole critical points structure looks like five 
straight lines running in all directions. An alternative 
technique to analyze this kind of problems is by using 



15 According to the implicit function theorem the equation y = f(x) 
can be inverted around the solution yg = f(xo) if /'(xo) 7^ 0. 
Critical points for which /'(xn) = define bifurcation branches 
if f"(x ) 5^0. 



the modular landscape of the exponential as explained in 
detail in Refs. [Hlls!- 

The upshot of this study is that e(u>) has only crit- 
ical points on the second Riemann sheet. This means 
that whenever the bound state becomes unbound we are 
crossing the branch cut, which happens for sufficiently 
small values of u. For the spectral integral that would 
set the condition M c < My /2 not to have cut intersection 
between the spectral function and the energy function. 16 



APPENDIX E: COMPUTING THE SPECTRAL 
INTEGRAL FOR THE VALENCE 
CONTRIBUTION FIRST 

In this appendix we show how our basic result for 
the valence contribution, Eq. (|86|) . can be reproduced 
by evaluating the spectral integral first. In fact, experi- 
ence with perturbative calculations of quantities such as 
/„■ or the pion form factor suggest that there is no danger 
in computing the ui integral first. At finite temperature 
this is actually the only way to get analytic results [88[. 
To compute the troublesome spectral integral first in the 
soliton case, we transform the scalar-pseudoscalar cou- 
pling into a vector-axial derivative coupling and a con- 
stant mass. We write 



1 - LOU 5 ) = U 5 



(El) 



where u 2 = U. This is a field dependent axial rotation 
and at the level of the determinant of the Dirac operator 
generates the anomalous Wess-Zumino term in the SU(3) 
case. Then, we get 



S aa >(x,x') = / dwp(u>) 



(E2) 



X (X 



u~ 5 (i$ + u^ lD i3u~ ti 



p(uj) 



u~ 5 \x') 



= [u( X y 5 (x\S (i@ + U- 5 i@u- 5 ) |a;Xa;')" 5 ] aQ; . 

Now, in the meson dominance model the propagator 
function is given by 

Sty) = [ doj^- = 2 mp {f) + i. (E3) 
Jc p-u p 

This relation is obtained for the closed contour under the 
assumption that p(u>) has a pole only at lo = 0, which 



16 Obviously, the smaller M c the better. Now, it would be interest- 
ing to see if there are chiral angles 9(r) for which one always has 
a bound state. In such a case cuts would not intersect. This is 
not such a strange situation, since e.g. in one dimensional quan- 
tum mechanics, any tiny attractive potential generates immedi- 
ately a bound state. Also, in 1+1 dimensions the chiral angle 
8(x) = asign(a;) has always a bound state e = Mcos(a) [87ll . 
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is clearly the case of the meson dominance model (1201) . 
Equation (|E3[) is just a functional relation, and does not 
depend on the fact that the variable j> is a matrix. We 
have indeed computed the spectral integral, but for this 
formula to be of any use, we have to look for the eigen- 
value problem of the chirally rotated Dirac operator 



u 5 it 



(E4) 



This eigenvalue problem may be tricky in practice since 
the operator is not normal, i.e. the adjoint operator does 
not commute with the original operator even if we go to 
the Euclidean space (where the derivative and the vector 
part of the coupling are anti self-adjoint and the axial 
part is self adjoint). This means that M n may be complex 
in general. 

Undoing field dependent chiral rotation we get 



where 



M n U 5 )y n = 0, 



(E5) 



(E6) 



Thus, the eigenvalue problem corresponds to the search 
for the (eigen)masses of the original Dirac operator which 
yield a zero mode of the four dimensional Dirac oper- 
ator. This equation looks as an on-shell condition for 
the bound quarks. In the stationary field case we have, 
^> n (x,t) — e~ iet \j} n {x). At this point we must make a 
choice on the values of e. A natural condition would be to 
request an exponential damping in the Euclidean space. 
Thus, we take e real and positive for forward propagat- 
ing states, t > 0, and negative for backward propagating 
states. (Another possibility that would be to take a con- 
tour with conditions on Ree), 



[-ia • V + M n (3U 5 ] ip n =eip n . 



(E7) 



Thus, here we encounter just the opposite problem to 
the case of the spectral mass problem; for a given en- 
ergy we have to look for the eigenmasses, M n (e) (instead 
of e n (uj)). Actually, this looks like an isospectral prob- 
lem; given the energy how many chiral couplings pos- 
sess the same energy. Thus, the eigenvectors of this 
problem are not exactly the eigenvectors of the standard 
Dirac Hamiltonian problem. Note also that because we 
are searching for zero modes of the Dirac operator its 
eigenvectors are also eigenvectors of the Dirac Hamilto- 
nian, which is not the case for non-zero eigenvalues, i.e. 
7o(«9t — H)^> n = X n ^n with A„ ^ cannot be solved by a 
eigenstate of H because of the 70 . Clearly, multivalued- 
ness issues may become quite relevant when discussing 
the possible equivalence of both problems. On the other 
hand, ffl = H if and only if M„(e) is real. For such a 
case if both M„ and e are real, then bound states hap- 
pen for — M n < e < M n . In general M n is compl ex, since 
even in the free case, U — 1, we have M„(e) = -y/e 2 — k 2 , 
with k n denoting the momentum quantized due to box 



boundary conditions. Then, we are lead to 
S aa >(x,x') = u- 5 (x)(x\S (i@ + u- 5 i$u- 5 ) \x') U - 5 {x') 

Hx) J2 M*)s ( M «) *J,0O«~ V) 



J2*n(x)S(M n )¥ n (x') (E8) 



+ e 



<*-*')(>(-, 



6{-e)6{t-t') S(M n (e))iPt(x'), 



In this representation, the states propagating forward in 
time have been chosen to be those of positive energy and 
we are led to the basic integral 



f°° de f° 
/ e~^S(M n {e))= / 

Jo 27r Jo 



* e -r £ +logS(M„(6)) / E9 % 

2tt 



with t — ► 00 . The function S(M) of Eq. (|E3| has no 
poles but branch cuts starting at M = zLMy/2. The 
stationary points are determined from the equation 

T= aio g S(M n (e)) =K{e) S>(M n M (E10) 



dc 



S(M n (e)) 



Note that a real S(M) requires M n < My/2, thus for 
M n (e) 7^ My/2 we must have a divergent M' n {e) — > ±00. 
Obviously, the smallest positive e (the valence orbit) ful- 
filling 1/M' n {e) — dominates, and we get the relation 



to = mm e 



M4(e)=oo, e>0 



(Ell) 



in agreement with our result of Sect. (JTTTJ) , which is the 
result equivalent to ([M]) written in terms of the inverse 
function. 



APPENDIX F: EVALUATION OF THE DIRAC 
SEA CONTRIBUTION TO SOLITON 
OBSERVABLES 

Dirac sea contributions to observables involve the com- 
plex mass integral. However, it is possible to rewrite this 
contribution as a real mass distribution. The vector spec- 
tral density py can be written as 



du 2 



where 



1 1 



2?ri 12w 5 (l-Aucu 2 ) 1 / 2 ' 



(Fl) 



(F2) 



The function py(w) has the property that its integrals 
with smooth functions vanish along the parts of the con- 
tour C encircling the branch points of Fig. [TJ as well as 
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are finite along the cut. Similarly, for the scalar part, ps, 
we have 



d 2 

Ps{u) = x^AsM 



with 



1 1 



u=l/Ml 



W 3 



27ril2w 4 Mj(l-4uw 2 ) 1 /2- 



(F4) 



We now define 



'(«) 



disc[py(w)] 



2 / dt^disc[ / o\/(u;)] ^ sca ^- — 
/i/(2V^) ' 2 

1 



12ttw 5 V4uw 2 - 1 ' 



(the factor of 2 comes from the two sections of C) and, 
similarly, 



A sea ,s(u) = 2 / duj&isc[ps{u, 
disc[p 5 (w)] - ^ 



(F6) 



The Dirac-sea contributions to observables are now ob- 
tained from 



(F7) 



«=1/Af? 



Thus, we have managed to "put the model on the real 
axis" , at the only expense of the need of differentiation 
with respect to u at the end of the calculation. As a bonus 
we also get a very high degree of convergence at large 
lo, since disc[py(w)] ~ 1/lu 6 and disc[ps(w)] ~ l/io 5 . 
It is useful to get rid of the integrable singularities in 
Eqs. (|F5IF6P by means of introducing the variable 



z = \J Aulo 2 — 1, 



2^ 



Then 



t , , 8u 2 
A S c a ,v{u) = — — 



dz 



1 



3^ Jo (z 2 + I) 3 

A SR 



(F8) 



(F9) 



2V^ 



and 

A S ea,s{u) 



N C My 



(4«) 



3/2 



o (^ 2 + 1) 5/2 



(F10) 



1 

x - 
2 



4, 



\/i T TT 

2^ 



-4 



2V^ 



The Dirac Hamiltonian has the form 

H(u>, to) = -ice ■ V + (3m + (3u)U 5 , 



(Fll) 



with the corresponding Dirac equation 

H(u>, m)\i; lo, m) — £i(uj, m)\i; lo, to) 



(F3) It has the property 



j 5 H(to,m)j 5 = H(-u,-m), 
j 5 \i;uj,m) = \i;-u>,-m). 



(F12) 



(F13) 



Thus 75 flips the sign of the current quark mass m and 
lo. Obviously, the spectrum of the operator is unchanged 
under this similarity transformation. Nevertheless, it is 
convenient in the numerical work to deal with positive 
values of to only. Thus, any one body observable at neg- 
ative lo can be transformed according to 



(F5) A sca (to) 



(i; to, m\a(u>, m)\i; lo, to) 



y^(»i m|757 5 a(cj, 771)7575 \i; lo, to) 



-w, -TO|7 5 a(w,TO)7 5 |i; -u, -to) 



(F14) 



Some examples of 750(0;, 771)75 encountered in the eval- 
uation of observables are 75^(0;, 771)75 = H(—lo, —to), 
75/375 = ~~ Pi e tc The formula for the energy following 
from the above formulation is given in Eq. (|106I107[) . 



APPENDIX G: INSTABILITY OF THE LINEAR 
MODEL 

Although SQM is only constructed in the nonlinear 
case, it looks tempting to extend it in the spirit of the 
original bosonizcd NJL model to a linear version where 
the fields may depart from the chiral circle, 

I = i dtop(Lo)Ti \og(i(^ — m — w(s + ij$T ■ p)) 
Jc 



^Jd^+p 2 ). 



(Gl) 



Here s and p denote the scalar-isoscalar and 
pseudoscalar-isovector fields, and g is a coupling 
constant. The meson fields s and p are dimensionless. 
In the chiral limit and in the vacuum (s = 1, p = 0) we 
find from the Euler-Lagrange equation for the s field the 
condition 



9 Jc W-l 



d 4 k 



-AiN c Nf diop(uj) j (27r)4fc2 _ (ws)2 . 



(G2) 



With the explicit form of the meson-dominance spectral 
function ([20|) and the techniques of Ref. [HH we find 



1_ = N c NfMy 

q 2 48tt 2 



(G3) 
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The effective potential assumes the form 



V 



fd*x 



N c NfMy 

1927T 2 



{2s 2 -s% 



(G4) 



where the s 4 term originates form the first term in 



Eq. (|Gip . This corresponds to an inverted Mexican Hat, 
and clearly displays instability. Therefore the linear ver- 
sion of the model (|G1[) does not make sense. 

It is "worthwhile to mention that this feature should 



not be regarded as specific to the spectral regularization. 
The inverted potential also arises when renormalizing the 
fermion determinant with the help of the ^-function reg- 
ularization in the local NJL model. Therefore one needs 
to use from the outset the non- linear realization of chiral 
symmetry on its own, and it cannot be treated as an ap- 
proximation to the linear model. In that regard we also 
note that the soliton instability in linear NJL models has 
been found in Ref. r 89L 
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